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Abstract 


Shock wave (Hugoniot), shock temperature, and release data are presented 
for several geophysically important, refractory materials. A sensitive multi- 
channel optical pyrometer was deveiopcvd to measure shock temperatures (2500 
to 5600®K at pressures from 48 to 117 G,?a) in anorthite (CaAlgSlgOa) glass. 
Shock temperatures of 3750 to 8000'’K at pressures from 140 to 182 GPa were 
measured in calcium oxide (CaO). Temperature data were used to constrain the 
energetics of the B1-B2 phase transition at 70 GPa in»‘CaO, and to construct a 
finite-strain equation of state for CaO consistent with previous Hugoniot data. 

The new CaO equation of state was used with equation of state parameters 
of other oxides to construct a theoretical mixed oxide Hugoniot of anorthite, 
which is in agreement with new Hugoniot data above about 50 GPa, determined 
using new experimental techniques developed in this study. The mixed oxide 
model, however, overestimates the shock temperatures, and does not accurately 
predict measured release paths. Both shock temperature and release data for 
anorthite indicate that several high pressure phase regions of stability exist 
above 50 GPa. A similar mixed oxide Hugoniot was constructed for lunar gab- 
broic anorthosite, and agrees with two new Hugoniot points at 120 GPa. Release 
data from lunar gabbroic anorthosite shocked to 120 GPa give evidence for 
shock vaporization. 

Because the densities and bulk properties of CaO and the high pressure 
phase or phases of anorthite are so close to those determined seismologically 


for the lower mantle, the amount of these materials present in the lower mantle 
is not well constrained. The possibility of significant enrichment of the lower 
mantle in these refractory materials, as predicted by inhomogeneous accretion 
models, is still open. 

A simple model is developed to explain the measured time dependences of 
radiated light in the shock temperature experiments, and constrain the absorp- 
tion coelTicient of the shocked material. The absorption coefficient is found to 
be an increasing function of shock pressure in shocked anorthite glass. 

Hugoniot and release paths were determined using electromagnetic parti- 
cle velocity gauges for San Gabriel anorthosite and San Marcos Gabbro shocked 
to peak stresses between 5 and 11 GPa. The data indicate a loss of shear 
strength in both rocks, and a partial phase transition of the anorthosite to a 
denser phase. This implies that estimates of shock wave attenuation in these 
materials based on elastic-plastic models are too high, and previously calculated 
amounts of internal energy gained by surface materials from impact or explo- 
sion events have been underestimated. 
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Introduction 


Dynamic compression of materials by strong shock waves provides a unique 
method by which matter can be studied In the laboratory at high pressures and 
temperatures. Shock waves have become a tool of both physicists and geophysi- 
cists, and furnish a research area rich in possibilities for experimental innova- 
tion. The purpose of this thesis is threefold: to develop new experimental shock 
wave methods to study condensed matter under high pressure, to explain the 
physics of experimental observations, and to apply the knowledge gained to 
problems in solid earth geophysics. 

In most experimental shock wave work, only the Kugoniot state is deter- 
mined. The Hugoniot curve, or locus of shock states of a material with a 
specified initial state, is useful in modeling high pressure bulk properties and 
detecting phase transformations with large volume changes, but its sensitivity to 
energies of transformation, specific heat, and Gruneisen parameter is limited. 
The usefulness of shock temperature measurements in constraining these ther- 
mal properties has been demonstrated, as well as the need for more precision 
and sensitivity in shock temperature measurement, A new optical pyrometer 
was designed with this purpose, and is described in Appendix I. Calibration, of 
this instrument is outlined in Appendix II, and determination of temperatures 
from measured spectral radiances, along with the data, is presented in Appendix 
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Chapter I illustrates the usefulness of shock temperature data in supple- 
menting Hugoniot data to constrain the physical properties of materials at high 
pressure. With both types of data taken together, the energy of the B1-D2 phase 
transformation of CaO is strongly constrained. Further constraint is provided 
for both heat capacity and Gruneisen parameter. Chapter 2 points out the appli- 
cation of temperature data to the detection of phase transformations which are 
unresolvable with Hugoniot data alone, in the case of anorthite. 

Experimental shock wave research at lOQ GPa pressures requires a large 
investment for little data. It is therefore important to glean as much informa- 
tion as possible from a given experiment. Chapter 3 contains an outline of three 
approaches used to obtain shock wave Hugoniot data from experiments with 
other primary objectives. 

Hugoniot data by itself is useful in determining high pressure bulk proper- 
ties of materials (i.e. bulk modulus, pressure derivative of the bulk modulus, and 
zero pressure density of high pressure phases). Chapter 1 provides a good 
example of an "inverse" problem— that is a least squares fit for these parameters 
to measured Hugoniot data with an assumed functional form for the high pres- 
sure behavior (third order finite strain). On the other hand. Chapter 3 illus- 
trates the "forward" problem, in which these parameters are calculated based 
on an a priori model— in this case a mixture of the oxide components of anorth- 
ite and lunar anorthosite— and a theoretical Hugoniot is cedculated. The validity 
of the model is checked by comparing the theoretical Hugoniot to the measured 
Hugoniot. 

Release data can supplement temperature data as indicators of phase 
transformations. Partial release states of shocked anorthite delineate possible 


regions of stability in Chapter 3, Partial release states of shocked lunar 
anorthosite indicate vaporization upon release. At lower stresses, release data 
can detect the loss of shear strength in shocked materials, as in Chapter 4 for 
San Gabriel anorthosite and San Klii?‘Cos gabbro, In Appendix IV, an experimen- 
tal technique is outlined which can determine an entire decompression path in a 
single experiment. This approach may eventually be applied to the problem of 
determining principal isentropes of fluids. 

The behavior of materials under high pressure is of geophysical interest for 
two principal reasons. First, the static pressures present in the earth's interior 
extend to 370 GPa, and knowledge of the high pressure physics of earth materi- 
als is a prerequisite to understanding the earth's interior. Second, because of 
the high dynamic pressures existing during the hypervelocity impact of bodies, 
an understanding of behavior of earth materials under shock is necessary in 
order to model impact and accretlonal processes. 

In Chapters 1 and 3, the bulk properties of CaO and of the high pressure 
phases of anorthite are compared to those determined seismologically for the 
lower mantle. Because of the similarity in bulk properties, the possibility exists 
for an enrichment of the lower mantle hi these materials, as predicted by inho- 
mogeneous accretion theories. 

The shock vaporization of lunar anorthosite suggested by the data in 
Chapter 3 has implications in cratering mechanics, and also supports the possi- 
bility of an impact-origin of the moon. The measurement of rarefaction wave 
velocities in Chapter 4 will allow more precise calculations to be made of shock 
wave decay in natural impacts, and more accurate estimates of energy parti- 
tioning in such events. 


Chapter 1 

Shock Temperatures in CaO 


ABSTRACT 

Blackbudy tp'^peratures of CaO shocked to pressures from 
140 to 182 GPa have been measured in the 3750 to 6000 ®K range 
using the Lawrence Livermore National Laboratory light gas gun, 

These shock temperatures, along with Hugoniot, porous Hugoniot, 
and isothermal data, are used to construct equations of state for 
the high pressure (B2) phase of CaO. The zero-pressure density of 
the B2 phase is between 3.8 and 4.0 Mg/m®, and the B1-B2 transi- 
tion energy is 2.1 to 2.3 kJ/g. The density and bulk modulus at 
pressures from 70 to 130 GPa are similar to seismically deter- 
mined values for the lower mantle. Thus, mixed-oxide models of 
the lower mantle will be insensitive to the quantity of CaO 
assumed, 

intiioduction 

The calcium-bearing minerals perovskite (CaTiOa), melilite 
(CaaAlgSiO^-CagMgSiaO^), diopside (CaMgSigOg) and anorthite (CaAl2Si20a) are 
among the first phases to condense from the solar nebula (Grossman and Lari- 
mer, 1974). They are of geophysical importance because they are hosts for the 



actinide series elemen*rtS U and Th, and because Ca Is the most common large- 
ion llthopbUe element in the solar system (Ross and Aller, 1976). In inhomogene- 
ous accreticm models (Turekian and Clark, 1969), the lower mantle Is 
significantly enriched in these minerals. At the high pressures of the lower man- 
tle (70 to 135 GPa), these minerals transform into a series of high pressure 
phases (Liu, 1976, 1979) whose high pressure equation of state properties may 
he close to those of their equivalent mixed oxides (Boslough at al„ 1983). Thus, 
the high pressure behavior of CaO is of interest. 

CaO exists In the B1 (NaCl) phase under standard conditions, and densities 
calculateu by extrapolation of its compression curve (e.g. Perez-Albuerne and 
Drickamer, 1985, Chang and Graham, 1977, Mammone et aL, 19B1) to lower man- 
tle pressures are significantly lower than seismologieally determined lower man- 
tle densities (e.g. Dziswonski at aL, 1P75, Anderson and Hart, 1976, Dziewonski 
and Anderson, 1981, Fig. 1<1). Jeanloz et al. (1979) found that at about 70 GPa, 
CaO undergoes a phase transition to the B2 (CsCl) structure, which is denser 
than the Bl phase. At high pressure the B2 phase of CaO has densities and 
compressibilities which are close to those of the mantle, according to Hugoniot 
data (Jeanloz and Ahrens, 1980). 

The reason for conducting shock temperature experiments on this material 
is to characterize its thermal properties at high pressures. Not only can the 
energetics of the Bl~B2 phase transition be constrained, but the possible 
existence of further phase transitions, such as melting, can be determined. The 
temperature-pressure Hugoniot is much more sensitive to such transitions than 

is the pressure-density Hugoniot (Lyzenga et of,., 1983). Armed with well- 
constrained transition energies in this way, the adiabatic Compression curve 
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Fig. t-1 Single crystal Hugoniot data (Jeanloz and Ahrens, 19B0), porous 
Hugoniot data (Marsh, I960), and static isothermal data (Mammone et al., 1981) 
for CaO. Solid and dashed curves are low and high zero-pressure density fits to 
B2 data, respectively. Calculated Bl isentrope based on ultrasonically deter- 
mined values of Kos= 112.5 GFa and K'og=4.8, of Chang and Graham (1977). Lower 
mantle curve based on PREM model of Dziewonslcl and Anderson (1981). 
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(Kos and K'os th® Birch-Murnaghan approximation, Davies, 1973) can be deter- 
mined with better precision, 

JEbq>erimeiaitai Efethods 

Tne CaO samples were obtained from Oak Ridge National Laboratory where 
they were grown by an arc-fusion process. They are from the same boule as 
samples used by Jeanloz and Ahrens (1980), and are free from all but trace 
impurities. Because optical methods are used to determine the shock tempera- 
ture, samples ideally should be perfectly transparent. These samples appear 
bluish in color, due to the Rayleigh scattering of white light by microbubbles 
inherently present in highly pure crystals grown by the arc-fusion method (Abra- 
ham et al., 1971). A Cary 17 spectrophotometer was used to measure the absor- 
bance at visible and near infrared wavelengths (Fig. 1-2). It demonstrates an 
approximate A* dependence, as is expected with Rayleigh scattering (Jackson, 
1975). 

There is some evidence for minute amounts of hydration in the form of 
portlandite (Ca(0H)2) in infrared spectra (Jeanloz and Ahrens, 1980). Both the 
microbubhies and the portlandite should lead to bulk densities which are lower 
than the x-ray density. This is observed: densities of samples used range from 
3.341 to 3.345 Mg/m^, compared to the x-ray density of 3.345 Mg/m®. Samples 
were approximately 3 mm thick and 17 mm wide, and both surfaces were pol- 
ished to optical quality. They were mounted on copper driver plates and 
shocked to pressures from 140 to 182 GPa by 2 mm thick copper and tantalum 
flyer plates fired from the Lawrence Livermore Laboratory light gas gun 
(Mitchell and Nellis, 1981) at speeds from 6.3 to 7.1 km/s. 


Absorbance 


Theoretical 
Rayleigh scattering 
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Measured CaO 


2bU 4UU b50 700 8 5 d ~IOOO 1150 I 30C 

f. , , Wavelength (nm) 

Fig. 1-3 ^.orbanc= spectrun. of a rapreseataUvo sampla of CaO. moasurod with 

h7- 1 =“rve is Of form: [h, . rfd with 

hgdo - 1.0, udo = 0.1, andX- = 400nm. A , 



The projectile velocity is determined by flash radiography, and this is used 
in an impedance match solution {Rice Bt eU., 1858) to calculate the shock pres- 
sure in the CaO. Thermal radiation is emitted from the shocked material and 
transmitted through the unshocked portion of the sample, which acts as a win- 
dow. The intensity of the li«ht is measured at six wavelengths in the ’asible and 
near infrared with an optical pyrometer (Lyzenga and Ahrens, 1979). The meas- 
ured spectral radiances are used to calculate the shock temperature. 

In an ideal situation in which the sample is initially perfectly transparent, 
the shocked sample is perfectly opaque, and th(»! shock wave is planar, parallel 
to the sample faces, and exhibits no time dependent behavior, the shock wave 
profile would appear as in Fig. 1-3. At time t^, the shock wave enters the trans- 
parent sample from the opaque driver plate. The shocked material immediately 
begins to radiate at an intensity determined by its temperature, and the signal 
increases suddenly to Rg. As the shock wave approaches the free surface, it is 
steady so the signal is constant. At time tg the shock reaches the free surface 
and reflects as a rarefaction. 'The light intensity drops to Rr, characteristic of 
the residual, post-shock temperature. 

Results and Discussioa 

A typical shot record from this study is shown in Fig, 1-4. It deviates from 
the ideal case in two principal ways. First, the rise from t^ to t^, is not instan- 
taneous, but is nearly 50 ns. Second, the light intensity is not constant during 
shock transit, but increases approximately linearly with time. 

The long rise time cannot be accounted for by the HP50B2-4207 photodiode 
rise time, which is approximately 1 ns (Lyzeno,® and Ahrens, 1979) nor is it due 


spectral Radiance 
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Time 

Fig. 1-3 Hypothetical shock temperature record for ideal case -where initial 
state is perfectly transparent and shock state is perfect black body. Shock 
wave enters sample at time tj and reaches free surface at Rs and are 
spectral radiances characteristic of shock temperature and post-shock tem- 
perature, respectively. 
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Fig. 1-4 Shock temperature record for shot Ca03T at X = 508nm. Time t,: shock 
wave enters sample, light is scattered indirectly to detector. Time t\,: shock 
wave enters detector field of view. Time t^: shock wave reaches free surface; 
spectral radiance is read at this point. 
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to shock wave tilt within the field of view of the photodiode. The active area of 
the detector is about 1 mm in diameter, and the optics demagnify the image by 
a factor of 0.43, so the field of view is 2.3 mm. The maximum projectile tilt is 
0.05 radians (Mitchell and Nellis, 1981). At typical projectile velocities of 
6mm/ jjLs, the maximum time interval for the line of intersection between shock 
front and driver/sample interface to cross the field of view is 20 ns; insufficient 
to account for the measured time lag. Because the samples are strong optical 
scatterers, the material within the detector field of view begins scattering light 
from the shock wave the instant it enters at one edge of the sample. The line of 
intersection between shock front and driver/sample interface can take as much 
as 70 ns to sweep across the 8 mm from the sample edge into the photodetector 
field of view, accounting for the long observed rise times. Unfortunately, 
interference from this scattered light precludes precise determination of the 
shock transit times, so these experim^^'nts cannot be used to determine 
Hugoniot states (e.g. Boslough et aL, 1983 for anorthite glass). 

The approximately linear increase in light intensity with time as the shock 
wave propagates through the sample can be explained by scattering and absorp- 
tion of light. The intensity of light transmitted through a medium is decreased 
by a factor of where is the linear absorption coefficient, h is the 

extinction coefficient due to scattering, and d is the light propagation distance 
through the medium. This distance is (t<.-t)Uj, where U, is the shock velocity, 
and tc is the time of shock arrival at the free surface. For (/z+h)d«l, the light 
intensity is 


I = Io(l - (m + h)(tc - t)U,) 


( 1 - 1 ) 


which increases linearly with time. At t5=to, the measured intensity is equal to 
the radiated intensity, so shot records were read at free surface arrival t,.. 


The ratio of initial intensity Ib to final intensity Ic was found to be smaller at 
shorter wavelengths, consistent with Rayleigh scattering. In Table 1-1, meas- 
ured values are tabulated for shot Ca03T, along with the best fitting theoretical 
values which assume a wavelengtfcu dependent Rayleigh scattering component 


h = ho 




and a wavelength independent absorption component (X. The effects 


of the two sources of attenuation are about equal at X=450 nm, and the 
wavelength independent component dominates at X=792 nm. 


Shock Temperaturea 

Two approaches were used to determine shock temperatures from optical 
pyrometer shot records. Both involve least squares fitting the spectral radiance 
data to a Planck distribution function (Appendix III). In method (a), both the 
emissivity and temperature are varied to obtain the best fit. In method (b) only 
the temperature is allowed to vary: the emissivity is assumed to be unity- 
equivalent to the assumption of a black body. 

Temperatures and emissivities calculated via methods (a) and (b) are 
presented in Table 1-2 with the standard deviations a and mean fractional devia- 
tions Pressures are calculated from an impedance match solution using the 

measured impact velocity and the known Hugoniots of the flyer materials 
(McQueen ef of., i970) and CaO (Jeanloz and Ahrens, 1980). In two of the three 
experiments the best emissivities are greater than unity— an unphysical situa- 
tion. In all experiments method (a) does not give significemtly better agreement 
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Table 1-1 

Shot Ca03T measured and theoretical intensity ratios 


X (nm) 

(maaaujrad 

^Ihacrstisal 

450 

.58 

.56 

508 

.46 

.46 

545 

.37 

.41 

598 

.42 

.37 

650 

.37 

.35 

792 

.29 

.32 


(mcBBUred ~ ^b/ 

r 1 

^thaoratlcal ~ |M 4 ho(Xo/ \)*jdj 


hodo=.27 

/zdg=.S9 


best fit to data 


Xg = 450 nm 
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Table 1-2 


Calcium Oxide Shock Temperature Data 


Shot 

Flyer /Driver 
Material 

Projectile 

Velocity 

(km/a) 

Initial 

Density 

(Mg/m») 

Presaure 

(GPa) 

Method 

Temperature Calculation 
c T a 

(“K) (W/srm3) 

a/\ 

Ca02T 

Cu/Cu 

6,l»e 

3.3449 

140.1 

(«) 

0.60 

4075 

5.6 

0.25 



;fc.06S 

±.0013 

±1,9 



±156 

X10,“ 







0>) 

'1.00 

3749 

7,8 

0.34 








±178 

X10.“ 


CaOlT 

Cu/Cu 

7.102 

3.3422 

171.0 

(«) 

1.09 

5242 

3.4 

0.21 



±.021 

±.0012 

±1.2 



±222 

X10.‘* 







(^) 

1.00 

5333 

3.5 

0.21 








i238 

xiO.*® 


CeOST 

Ta/Cu 

6,457 

3.3411 

181.6 

(«) 

1.15 

5699 

4.2 

0,14 



±.027 

±,0014 

±1.2 



±195 

X10.‘* 







(b) 

1.00 

6083 

4.5 

0.15 








±225 

xio.« 



(a) : best fit with Planck function and variable s 

(b) : best fit with Planck function and e=l (blackbody) 



than method (b), as measured by the standard deviation. Shot Ca02T was the 
only experiment which led to a best fit with £<1. If the emissivity were really 
0,6, as determined by the flt, one would expect to see the time dependent 
behavior similar to that observed in forsterite by Ahrens at al. (1982) or anorth- 
Ite glass at lower pressures (Boslough et al., 1903), due to the increasing thick- 
ness of the shocked material with time, unless the shock front is opaque and has 
a finite reflectivity. Because this is not observed, it is reasonable to conclude 
that the emissivity is close to unity for all three experiments. 

Equation of State 

Shock temperatures determined by methods (a) and (b) are plotted in the 
P-T plane in Fig. 1-5, along with a theoretical calculation by the method of 
Ahrens et al. (1969), based on parameters of Jeanloz and Ahrens (i960), and a 
Dulong-Petit specific heat Cy= 3R. If the temperatures calculated by method (a) 
are believed, there is strong evidence for a phase transition at about 4000‘‘K, 
with a transition energy, Eir&0.5 kJ/ g. We do not favor this interpretation, but 
include it for compl'?.teness. 

The results of method (b) (blackbody temperatures) are preferred for rea- 
sons outlined above, and there is no evidence for a phase transition at or above 
4000"K based on these temperatures. However, all three measured tempera- 
tures lie about 600®K below the theoretical curve of Jeanloz and Ahrens, Indicat- 
ing that Etr, the energy of transition between the initial state and the metastable 
final phase at standard conditions, is actually larger than the value used in the 
calculations: 1.6 kJ/g. 


Temperature (“K) 


6000 


^000 


o method (a) 




150 

Pressure (GPa) 


« 


X 


200 







- 18- 


f>p 


One possible Interpretation Is a second phase transition below 140 GPa. 
This phase transition would have to have an Eir«0i5 kJ/gr unless Ejf tur the B1 to 
B2 transition Is significantly less than 1,6 kJ/g, which is unlikely. Such a small 
Eir would probably not result from melting, which would be expected to have an 
Etr of around 1,4 kJ/g, the latent heat of fusion at 1 atmosphere (Robie et al., 
1979). 


If the B2 phase indeed exists on the Hugoniot to 102 GPa, a complete equa- 
tion of state for that phase can be constructed from the present data. Rather 
than simply assuming the parameters determined by Jeanloz and Ahrens and 
finding the Eu- which provided the best agreement to the new temperature data, 
all the available data for the H2 phase were simultaneously fit. This included 
porous Hugoniot data (Marsh. 1980), and Isothermal data (Mammone et of., 
1901). The inversion procedure was based on third order finite strain theory 


(Davies, 1973) and the Mie-Gruneisen approximation, with ^assumed constant. 


where y is the thermodynamic Gruneisen parameter V 


3P 


^as jv 


. Despite the abun- 


dance of data, a wide range of values provided good conformity. Parameters for 
the two extreme models are given in Table 1-3, based on zero pressure densities 
of 3.8 and 4.0 g/cm®. 

Curves calculated with these parameters are plotted with the shock tem- 
perature data in Fig. 1-5, and with the isothermal and Hugoniot data in Fig. 1-1. 
Both the low and high density solutions give approximately the same curve in 
the P-T plane, but the high density equation of state requires a 0.2 kJ/g greater 
' energy of transition. The slope of the theoretical curve is determined primarily 
by the specific heat. It is in good agreement with the data, indicating that the 
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Dulong-Petit value (C^ =3R) is a good approximation for CaO, in contrast to 
anorthite (Boslough et al„ 1983) and quartz (Lyzenga et al„ 1963). 

Conclundoni 

The first shock tenciperature data have been presented for CaO to 182 GPa. 
They can be fit well with theoretical shock temperature calculations for CaO in 
the B2 phase and allow a precise energy of transition (2.1 to 2.3 kJ/g at STP) to 
be determined for the B1-B2 phase transformation. This Et, in turn provides 
further constraint for calculating the bulk modulus and its pressure derivative 
, from Hugoniot data. A range of equations of state were determined in this 
manner. The equation of state with the lowest zero pressure density in the 
range does not differ greatly from that of Jeanloz and Ahrens (1980) but is 
better constredned. At lower mantle pressures, the present models predict den- 
sities of CaO to be lower than the lower mantle by up to 1.3%, and bulk moduli 
greater by 3 to 13%. Thus, depending on the bulk properties of the other com- 
ponents of lower mantle models, and the particular model used, CaO cannot be 
precluded as a major component of the lower mantle on the basis of its densitj*- 
and bulk modulus. However, because such a wide range of equations of state for 
the B2 phase are indistinguishable with the data, theoretical considerations may 
be necessary to choose the best equation of state. 
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Chapter 2 

Shock Temperatures in Anorthite Glass 

ABSTRACT 

Temperatures of CaAlgSigOa (anorthite glass) shocked to pres- 
sures between 48 and 117 GPa have been measured in the range 
from 2500 to 5800®K, using optical pyrometry techniques. The 
pressure dependence of the shock temperatures deviates 
signiflcanlly from predictions based on a single high-pressure 
phase. At least three phase transitions, at pressures of about hh, 
85, and 100 GPa and with transition energies of about 0.5 kJ/g each 
(wi.5 kJ/g total) are reqiaired to explain the shock-temperature 
data. The phase transition at 100 GPa can possibly be identified 
with the stishovite melting transition. Theoretical models of the 
time dependence of the thermal radiation from the shocked 
anorthite based on the geometry of the experiment and the 
absorptive properties of the shocked material give good agreement 
with observations, indicating that it is not necessary to invoke 
intrinsic time dependences Lo explain the data In many cases. 
Observed time dependences were used to calculate absorption 
coefficient of the shocked material of from 1 to >14 mm~^; an 
increasing function of shock pressure. The assumption that the 
shocked material radiates as a blackbody is supported by the 


theoretical model, and by the close agreement between measured 

and calculated blackbody spectral radiance as a function of 
wavelength. 

Iniroductioa 

Because of the cosmochemical abundance of Ca and Al, and the refractory 
nature of calcium-aluminum silicates such as anorthite, their high pressure 
polymorphs are important to consider in any model of the earth’s mantle. 
Models for the earth’s composition based on chondritic abundances predict 

atomic abundance ratios of m 0.05 and 0.07 (Ross Sc Aller, 1978). It is 

Si Si 

possible however, that the lower mantle is significantly enriched relative to the 
whole earth in calcium-aluminum silicates, because these are among the first 
phases to condense from the solar nebula (Grossman & Larimer, 1974) and 
accrete onto the earth, in inhomogeneous accretion models (Turekian Sc Clark, 
1969). Moreover, the construction of equations of state for anorthosite are key 
ingredients in carrying out cratering calculations for impacts on the crusts of 
terrestrial planets (Ahrens Sc O’Keefe, 1977, 1983), and in determining quantities 
of melting and vaporization resulting from such an impact (Boslough Sc Ahrens, 
1903). 

In previous work on anorthite, a high pressure equation of state to pres- 
sures greater than 120 GPa has been constructed by JeanLoz Sc Ahrens (1980) 
based on Hugoniot data only. This equation of state is incomplete in that the 
thermal behavior of the shocked material was inferred from pressure-density 
states only. The parameter which is most sensitive to thermal properties and 
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energetics of phase transitions, namely the temperature, was not measured. 

Past work (Kormer et al„ 1965, Ahrens ef oZ., 1982, Lyzenga et al,, 1983, 
Bosiough et ol,, 1983a) has demonstrated the usefulness of shock-temperature 
data in constraining the high-pressure-high-temperature properties of materi- 
als, and in defining phase transitions. For example, the stishovite melting curve 
■was determined from shock-temperature data on a-quartz and fused quartz. 
This phase transition is readily apparent from the temperature data, but is not 
at all obvious from Hugoniot data in the pressure-density plane (Lyzenga et al., 
1983). It is reasonable to expect similar results for other materials. The pur- 
pose of the present study was to measure temperatures of shocked anorthite 
glass using the optical pyrometry technique, and determine the high-pressure 
thermal behavior of CaAlgSigOg. 


£q;)enmei]ital 

Amorphous CaAl^SigOg samples, with measured Archimedean initial densi- 
ties of 2.69 Mg/m^ were used in these shock-temperature experiments. The 
samples, obtained from Corning Glass Co. (Corning, N.Y.), were clear and homo- 
geneous, with a composition described by Angg g Abo ,4 Org i — essentially the com- 
position of pure anorthite, A microprobe analysis of a representative sample is 
given in terms of oxides in Table 2-1. The transmittance of a 3mm thick, pol- 
ished sample of the glass over the wavelength range of interest is shown in Fig. 
2-1, measured with a Cary Model 17 spectrophotometer. 

Shock experiments were carried out on two-stage light-gas guns at the Cali- 
fornia Institute of Technology and Lawrence Livermore National Laboratory. 
Copper and tantalum projectiles were accelerated to velocities between 3.7 and 


Oxide 


Wt.% 


NagO 

0.05 

MgO 

0.00 

AlgOs 

36.36 

SiOg 

43.48 

KgO 

0.01 

CaO 

20.73 

TiOg 

0.06 

FeO 

0.04 

BaO 

0.03 


Total 100.76 
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Fig. 2-1. Transmittance spectrum of a representative sample of CaAlgSigOa glass 
(sample An7T), measured by Cary 17 spectrophotometer. Transmittance shows 
negligible wavelength dependence between wavelength limits of optical pyrome- 
ters (dashed lines). 
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5.7 km/s. Projectile velocities were determined using flash X-radiography, as 
described elsewhere (Jeanloz & Ahrens, 1977; Mitchell &: Nellis, 1981). Shock 
pressures were determined by the method of impedance matching (Rice ei al., 
1958) and using the standard equation of state of McQueen et al. (1970) for tan- 
talum and copper, and the Hugoniot of Boslough et al. (l9B3b), for the anorthite 
glass. 

Temperatures were determined by measuring the intensity of thermally 
radiated light at several wavelengths in the visible and near infrared 
wavelengths. The thermal radiation was filtered by interference filters and 
focused onto silicon photodiodes. The resulting electrical signals were recorded 
by an array of oscilloscopes with Polaroid cameras. The recorded time-resolved 
voltages are directly proportional to the spectral radiances of the thermal radi- 
ation. Calibration of the pyrometer systems waa carried out using light sources 
with known spectral radiances in the wavelength range in which measurements 
were made. The two optical pyrometers used in this study are described in 
detail by Lyzenga & Ahrens (1979) and in Appendix 1 and II. 

Results 

Eight shock temperature experiments were conducted on anorthite glass; 
three using a 10 nm bandwidth six-channel optical pyrometer at Lawrence Liver- 
more National Laboratory (Lyzenga ic Ahrens, 1979), and five using a 40 nm 
bandwidth four-channel optical pyrometer at the California Institute of Technol- 
ogy (Appendix I). Representative data are shown in Figs. 2-2 and 2-3. 

Three general types of time-dependence were observed. In shots AnBT, 
An5T, and An3T~the three experiments below 70 GPa— the sharp increase in 


200ns 


100ns 


OF POOil QIJ 


200mV 


Fig. 2-2 Shock temperature records for CaAlgSigOfl glass. Signal voltage is 
directly proportional to the intensity of the radiated light at the wavelength of 
the Alter, (a) An5T; P = 55 GPa, T = 2700*K, X = 750nm. Time dependence is 
due to combination of ejects: Cooling of interface and thickening of radiating 
(shocked) layer, (b) AnlOT: P = 117 GPa. T - 5600“K, X = 800nm. Initially con- 
stant signal due to high opacity of shocked layer. Decrease in intensity just 
before free surface arrival due to decay of shock by eige rarefaction. 









Signal Voltage 


Of- 



)mV 


Time 

Fig. 2-3. Shock temperature record for shot An9T: P=85 GPa. T=4000“K, 

X=75Unm. Time dependence can only be explained by unsteady shock. Signal 
voltage is directly proportional to the intensity of the radiated light at 750nm. 
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intensity corresponding to the entrance of the shock wave into the sample is fol- 
lowed by an exponential decay to an approximately constant value (Fig, 2-2a). 
In shots An2T, An7T, AnlT, and AuiOT~four experiments above 09 GPa~the light 
intensity remains constant to within about 10% until the shock wave reaches the 
free surface or is" overtaken by an edge rarefaction (Fig. 2-2b), Shot An9T-to 
about 05 GPa— shows an anomalous behavior, in which the intensity decreases at 
a nonuniform rate during shot transit (Fig. 2-3). 

There are three causes which could account for time dependence in the 
emitted light; temperature dependence, emissivity dependences, or time 
dependence of the absorption or scattering of an intervening layer. All three of 
these causes are invoked to explain the observed time dependences. It is impor- 
tant to note that the observed light Is not necessarily radiated from only the 
shock front unless the shocked material is perfectly opaque. In general, a 
material that is irdtially semi-transparent, with Initial linear absorption 
coefficient a^ is shocKred to a state with linear absorption coefficient a, and tem- 
perature T. The absorbance of the unshocked layer is Ay = e"****^, and the 
shocked layer is Ag = e where d„ and d, are the thicknesses of unshocked 
and shocked layers, respectively. According to Kirchoff s radiation law for a 
non-reflector (Jenkins & White, 1976) the emissivity of a layer is equal to its 
absorbance. If the sample thickness is d, the shock velocity is U„ and graybody 
light (with a Planckian wavelength distribution and an emissivity less than unity) 
is radiated only from the shocked layer at constant temperature T, the time 
dependence will be 

I.(t) =fA(T) 


( 2 - 1 ) 


where !>(?) Is the Planck distribution function. In the case of anorthite glass, 
the initial state is essentially transparent (Fig. 2-1), so a„=0 and this becomes: 

= (2-S) 

It is also necessary to coaisider the effect of light generated at the interface 
between the driver and sample material, due to the non-ideality of the surfaces, 
which, despite their opticaj, quality, have imperfections which leave a gap on the 
order of l/xm wide as determined by counting visible fringes (Newton rings) with 
an optical flat. 

The temperature at the interface due to passage of the shock wave is higher 
than the shock temperature of either boundary material, due to multiple rever- 
beration (Urtiew ic Grover, 1P74) and shock heating of any gas trapped in the 
gap. The temperature decay as a function of time at the interface can be 
modeled simply as a Fourier diffusion boundary value problem in one dimension. 
The metallic driver is considered to be a thermally conducting half-space with 
initial temperature Td. the shock temperature of the driver, and thermal 
diffusivity ic. The sample is modelled as a thermal insulator. This problem is 
solved by Carslaw & Jaeger (1959, p.5l). The solution is 

where t is the time after impact, x is the distance of a point in the driver from 
the interface, Q is the quantity of heat per unit area (unknown) deposited at the 
interface by the closing of the gap, and p and Cy are the density and the specific 
heat at constant volume of the driver material, respectively. The time 
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dependence of the interface temperature is 

T,(t) = T(0,t) = Ta + Ct-^^2. (2-4) 

If the interface radiates as a blackbody’, the observed light intensity from it 
(with attenuation due to the growing shocked layer) is 

I|(t) = fx(Td + ( 2 . 5 ) 

The total light intensity observed is due to both sources (Eqns. 2-1 and 2-5) 

I(t) = f^(T)[l-e"“*^*^) fA(Ta+Ct-‘^ 2 )e“‘*^*\ ( 2 - 8 ) 

One can easily see that for t » — 77 -, the shocked layer is effectively an 

a,Ug 

opaque blackbody, and there is no time dependence in the observed intensity. 
For sufficiently large a,, this occurs at t smaller than the rise time of the pyrom- 
eter photodiode and no time dependent behavior is observed at all. This is 
observed in shots AnlT, An7T, and AnlOT (Fig. 2 - 2 b). 

For smaller a,, the time dependence is depicted schematically in Fig. 2-4, 

and is illustrated by the spectral radiance observed in shot An3T in Fig. 2-5. The 

contribution due to the shocked sample increases from zero, when the shock 
wave enters the s^ple, asymptotically to the intensity of a blackbody at 

temperature T, with a time constant r= - 7 , ' ■ The interface t ontribution 

decreases from infinity (because of the idealized assumption that all the heat is 
deposited at the two-dimensional interface) to zero, with a time constant r < 

— because the source is also decaying with time. The observed light 

a,U, 
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Fig. 2-4 Schematic of shock temperature experiment. Measured light intensity 
is sum of light from interface (Ij) and shocked sample layer (Ij. (a) Just after 
shock enters sample, light from hot interface dominates, and is not attenuated 
significantly by thin shocked layer, (b) Ij has decreased due to cooling of inter- 
face and thickening of attenuating layer. Ig has increased due to thickening of 
emitting layer, (c) Interface has cooled and is further blocked by thick 
shocked layer. 1,. the light intensity characteristic of the shock temperature 
now dominates. 
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CaAIgSigOg Glass 
Shot An 3T 650 nnn 


Observed 


'Shocked sample layer 
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Time (ns) 

Fig. 2-5. Measured spectral radiance from shot An3T (P = 66 GPa, 
T = 3000"K) at X = 650nra. The measured time dependence can be fit by a func- 
tion of the form: A(l— e~®^) + e an approximation of equation (2-6), 

with A = 4.3x10*^ VTSr~*m~^ B = 8.2x10*^ WSr-^m"*^, and a ~ 8.0xl0“^ns~^ 
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intensity is the sum which decreases asymptotically to f;^(T), with an effective 


decay time bounded above by 


a.U, 


This iu the time dependence observed in 


shots An2T, An3T, An5T, and An0T (Fig. 2-2a). 

One experiment, An9T, exhibits time-dependent behavior not accounted for 
by this sii-nple model. The intensity of the emitted light decreases in steps as 
the shock wave passes through the sample, indicating that the time-dependence 
is in the shock front itself, as a fluctuation in temperature, absorption 
coefficient, or both. This is similar to behavior observed in fused quartz shocked 
to the stishovite-liquid mixed phase region by Lyzenga et al. (19B3), which they 
interpret as a possible rate dependent behavior due to metastability near the 
phase transition pressure. 

In the simple model developed above (Eq. 2-6), the linear absorption 
coefficient a, for the shocked anorthite can be constrained by measuring the 
effective decay times from the oscilloscope records for the shots whose results 

1 


agree with the model. The upper bound for a, is a. 


smaz 


-. These values are 


tabulated in Table 2-2. The bound on a, is seen to be an increasing function of 
shock pressure. 


CaAIgSigOe i%ock Temperatures 

Temperatures were calculated from the spectral radiance data by fitting 
them to a Planck function, as outlined in Appendix III. This was done in two ways 
for each experiment. First, a least squares fit is carried out in which both the 
temperature T and emissivity e are independent variables. The temperatiares 
and emissivities determined in this manner are given in Table 2-3, along with the 


Table 2-2 


Bounds on absorptivity of shocked anorthite glass 


Shot 

T 

(ms) 

u. 

(mm//us) 

^imu 

Pressure 

(GPa) 

An6T 

0.15 

6.5 

1.0 

50.3 

An5T 

0.10 

7.0 

1.4 

57.0 

An3T 

0.08 

7.7 

1.6 

68.2 

An2T 

0.01 

8.8 

12 

90.4 

Others 

<0.007 

<10. 

>14 

>93.1 
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Table 2-3 

Anorthite Glass Shock Temperature Data 


Shot 

Flyer/ 

Driver 

Material 

Projectile 

Velocity 

(km/sec) 

Initial 

Density 

(H«/m3) 

Pressure 

(GPa) 

Temperature Calculation 
Method c T a 

(®k) (w/arm^) 

An6T 

Cu/Cu 

3.72 

2.6908 

48.3 

(a) 

1.76 

2323 

0.2 

0.02 



i.lO 

±.0010 

±2.2 



±5 

klO.io 







(b) 

1.00 

2481 

1.4 

0.13 








±37 

xl0.“> 


An5T 

Cu/Cu 

4.03B 

2.6909 

55.0 

(a) 

0.96 

2723 

2.7 

0.12 



±.002 

±.0013 

±1.1 



±37 

xlO.io 







(b) 

1.00 

2714 

2,8 

0.12 








±39 

xlO.‘o 


An3T 

Cu/Cu 

4.536 

2.6900 

86.4 

(a) 

1.05 

2997 

1.5 

0.30 



±.006 

±.0010 

±1.2 



±104 

xlO.ii 







(b) 

1,00 

3017 

1.5 

0.30! 








±105 

xl0.i» 


An9T 

Cu/Cu 

5.255 

2.6953 

84.5 

(a) 

1.01 

4002 

1.2 

0.03 



±,003 

±.0012 

±1.3 



±21 

xl0.i‘ 







(b) 

LOO 

4010 

1.3 

0.04 








±22 

xl0.“ 


An2T 

Cu/Cu 

5.455 

2.6899 

89.9 

(a) 

1.34 

3802 

0.4 

0.01 



±.010 

±.0011 

±1.3 



±6 

xlO.ii 







(b) 

1.00 

3985 

4.0 

0.12 








±71 

xl0.“ 


AnTT 

Cu/Cu 

5.563 

2.6927 

92.9 

(c) 

1.00 

4564 





±.006 

±.0011 

±1.2 



±104 



AnlT 

Ta/Ta 

5.678 

2.6893 

111.5 

(a) 

1.24 

4904 

1.1 

0.08 



±.010 

±.0011 

±1.3 



±76 

Xl0.l2 



• 




(b) 

LOO 

5ii0 

1.6 

0.12 








±122 

XIO.12 


AnlOT 

Ta/Cu 

5.562 

2.6944 

116.8 

(a) 

0.91 

5741 

3.5 

0.19 



±.003 

±,0012 

±1.1 



±218 

xl0.« 







(b) 

1.00 

5625 

3.6 

0.20 








±214 

xl0.»8 



(a) : best fit with Planck function and variable s 

(b) : best fit with Planck function and esl (blackbody) 

(c) : brightness temperature 





standard deviation a and the mean fractional deviation of the fit. There 

appears to be no correlation between the shock pressure and the emissivity 
which gives the best fit, and in the majority of cases the best fitting emissivity is 
greater than unity--an unphysical situation. 

Because of this discrepancy, and the fact that the model developed in the 
previous section implies that the observed intensity corresponds to a blackbody 
on the asymptote, temperatures were also calculated assuming e = 1.0, and con- 
ducting a standard single variable least squares fit for T. 

In four of the eight anorthite shots, the best fitting emissivities are within 
10% of unity, and the goodness of the fits, as determined by the standard devia- 
tions are not significantly reduced by requiring that s = 1.0. In three shots, the 
best fitting emissivities are significantly greater than unity. Two of these shots, 
An2T and An6T, have extremely good fits, with mean fractional standard devia- 
tions of .01 and .02, respectively. However, the emissivities which give rise to 
these fits are 1.34 and 1,76, respectively; unphysical values which imply that the 
quality of the fits is fortuitous. When the inversion is carried out with an emis- 
sivity of unity for these three shots, however, the qualities of the fits are still 
comparable to thoj;e for the other experiments. 

In one shot, (An7T), the spectral radiance was obtained at only one 
wavelength. The shock temperature of this shot was assumed to be the bright- 
ness temperatuure. This assumption is based on the close agreement between 
the observed spectral radiances and blackbody distributions in the other seven 
experiments. The temperature error for this shot was taken from the estimated 
error in the measured spectral radiance. 


In the other seven shots, the temperature errors were determined sys- 


tematically by the equation 


” fT (T.\ = GOOnra) ’ 

where £7^ is the standard deviation in spectral radiance, and f>r(T,X) the partial 
derivative 'of the Planck distribution function with respect to temperature. The 
partial derivative was evaluated at a wavelength of 600 nm, because this is near 
the middle of the widest wavelength range used. 

The intensities used in the temperature calculations were taken from the 
oscilloscope records as indicated in Figs. 2-2 and 2-3. According to the simple 
model, the intensity of light most representative of thermal radiation from the 
shocked region of the sample is just before the shock wave arrives at the free 
surface. In all but one case, that is where the oscillograms were read. In the 
other case, shot An9T the shock wave begins to decay before it approaches the 
free surface, apparently due to an edge rarefaction. Thr? intensity was taken 
from this set of shot records just prior to the beginning of the decay of the 
shock wave. The spectral radiances determined in this manner are tabulated in 
Appendix III, along with plots which show the measured spectral radiances and 
those calculated for the blackbody temperatures which give the best fits. The 
results of shot An9T are illustrated in Fig. 2-8. 

Temperatures determined by both of these methods, along with their 
respective emlssivities, standard deviations in spectral radiance and fractional 
standard deviations are summarized in Table 2-3. The associated pressures are 
determined from the measured projectile velocity, and impedance matching 
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yrith a fit to the U^-Up data obtained in the Hugoniot experiments on anorthite 
glass (Boslough et ai„ 1983), These impedance-match determined pressures are 
the preferred values to use when modelling the pressure-temperature Hugoniot, 
despite the fact that pressures were determined more directly for each shot by 
measuring shock-transit times. In the more direct pressure determination, pro- 
jectile and shock velocity must both be measured. The precision by which the 
shock velocity can be measured is not as high as that of the projectile velocity 
measurement, and larger random errors can be introduced. In the impedance- 
match method, only the projectile velocity is used, which can be measured to 
better than ±0.2%, and a curve which is determined b 5 »^ a least-squares fit to the 
direct pressure measurements. This in effect smooths out the random errors. 

Diacussioii 

The measured shock temperatures in anorthite glass are plotted as a func- 
tion of Hugoniot pressure in Fig. 2-7, along with pressure-tempsrature Hugoniot 
calculations based on several models. These models are the Jeanloz Sc Ahrens 
(1980) equation of state with a Dulong-Petit specific heat (Cy = 3R), the same 
equation of state with variable specific heat, and a model in which several phase 
tremsitions occur. The latter two models can pro^de reasonable fits to the data. 
The method of Ahrens et al. (1969) was used to calculate shock temperatures in 
all cases. 

The Jeanloz & Ahrens (1980) equation of state considers all states above 30 
GPa on the Hugoniot of anorthite crystal to be a single high pressure phase, as 
found by McQueen et al. (1967), and fits the Hugoniot data with a third order 
Eulerian finite strain (Burch-Murnaghan) isentrope (Davies, 1973). The 


ORlGSf^ML PASlS 

OF POOR QUALITY 



Fig. 2>7, Measured shock temperatures of CaAlgSigOa assuming t - 1.0, and four 
calculated pressure-temperature Hugoniots. (A) Mixed oxide model of Boslough 
et al. (1983). (B) Jeanloz and Ahrens (1980). (C) This study, assuming, three 
phase transitions to high-pressure phase regimes denoted by I, II, and III, and 
Cy = 3R. (D) This study, assuming no phase transitions and Cy = A + BT. The 
lower two dashed lines are hypothetical phase transitions. The upper one is the 
stishovite melting curve (Lyzenga et al., 1983). Three shock temperature data 
by Schmitt and Ahrens (1983) are included, which exhibit heterogeneous hot 
spot (graybody) behavior. 
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Grunelsen parameter calculated by the pressure offset of the porous 
Hugoniot. The transition energy between single crystal anorthite and the meta- 
slable high pressure phase at standard conditions (Etj.) hi this equation of state 
is .72 MJ/kg. This corresponds to an E^j. of .46 MJ/kg for anorthite glass, with the 
.26 MJ/kg energy difference between anorthite glass and crystal is accounted for 
(Robie et al., 1978). Jeanloz & Ahrens considered anorthite above 2000°K to be a 
Debye solid in the classical limit— anorthite having a Debye temperature of 
1000'’K— and used the Dulong-Petit specific heat. This calculation gives reason- 
able agreement to the measured points below 80 GPa, but overestimates the 
shock temperature at higher pressure by as much as 1000®K Clearly, the slope 
described by the data is less steep than the calculated slope. 

In calculating shock temperatures by the method of Ahrens et al. (1969), 
the only way the theoretical T-P slope can change is to have a pressure or 
temperature-dependent Cy and/or phase transitions. There is ample theoretical 
and experimental evidence for a temperature dependent Cy in this temperature 
range, in departure from the Debye model of solids. A specific heat of the form 

Cy = A + BT (2-Q) 

is equivalent to that of Wallace (1972) for NaCl, where A=s3R and fl = Pg - 2Ag. 
Here, Pg is the electronic contribution and Ag is the first order anharmonic lat- 
tice contribution. The same form has been shown by Lyzenga et al. (1983) to 
describe shock temperature dependence on pressure for solid and liquid stisho- 
vite. 

The coefficients for the best fit to the anorthite data are A = 0.92 kJ/kg°K 
and B = 9.3x10“® kJ/kg®K®. These coefficients are comparable to those 
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determined by Lyzenga at al. (19B3) for solid and liquid stishovite, which has 
approximately the same mean atomic weight per atom as anorthite, with mean 
values of A = kJ/kg“K and B = 1.17x10 * kJ/kg®K^. B is In fact smaller for all 
the anorthite data taken together than for the single phases of the S 1 O 2 data, 
which indicates less temperature dependence for Cy overall. 

H is probably unreasonable to expect that the anorthite undergoes no 
phase transition along this region of its Hugoniot, however, when one considers 
the complexity of its phase diagram at lower pressures, below 30 GPa (Liu, 
1978), and when the behavior of simpler silicates such as SiOj in the same pres- 
sure range (Lyzenga at al., 1903) is considered. Anorthite disproporllonates to 
grossular, kyanite, and quartz at pressures from 3 to 15 GPa (Goldsmith, 1980; 
Boyd ie England, 1961) and at 15 GPa the quartz transforms to stishovite along 
its Hugoniot (McQueen at al. 1963). Phase changes in grossular have been 
observed above 25 GPa (Liu, 1979). Hugoniot data indicate a mixed phase region 
in anorthite form 10 to 33 GPa (Jeanloz k Ahrens, 1960). Jeanloz 8c Ahrens give 
four candidate high-pressure phases above 33 GPa: CaAlgSigOg (hollandite), CaO 
(B2 structure) + AlgOg (corundum) + 2Si0a (stishovite)— (mixed oxide assem- 
blage), CaAl 204 (calcium ferrite structure) + SiOg (stishovite), and CagAlgSiO^ 
(sodium titanate structure) + AlgOs (corundum) + SSiOg (stishovite). It is not 
unlikely, however, that phase transitions occur between these candidate phases, 
or others, above 33 GPa, or that further disproportionations or polymorphic 
transitions occur in the components of the candidate assemblages. Melting is 
also expected to occur in this shock-pressure range. 

With each phase change is associated an energy of transition, which can be 
expressed in terms of AEtri the energy diflerence between the metastable phases 


at standard conditions. The measured temperature for anorthite shocked Into a 
mixed phase region will, if In equilibrium, lie on the coexistence curve. For 
anorlhite shucked into a single phase regime, the measured temperature will lie 


^tr 

Cy 


below th«i;s calculated metastable temperature of the lower phases. 


The 


slopes of the single>phase segments of the Hugonlot depend on Cy, 


For constant Cy = 3R, at least three phase transitions (to assemblages indi- 
cated as I, II, and III in Fig. 2-7) are required to give a good fit to the data. If 
stishovite is a component of assemblage II, the highest transition~at about 

9 

4800°K~may correspond to the stishovite melting coexistence curve determined 
by Lyzenga et al. (1903) who found AEt, for the stishovite to liquid transition to 
be about 1.6 MJ/kg. The transition at 4800*^ in anorthite gives a temperature 
difference of about 450“ between the phases, which leads to AEt,^.5KCJ/ kg. The 
mass fraction of SiOg in CaAlgSi^Og is about .38, so if this pbasr^. iransition is melt- 
ing of SiOg only, this AEtr corresponds to 1.4 MJ/kg in the SiOg, which compares 
well with the value of Lyzenga at al. (1983). 


Conrliisions 

When shock-temperature data are included in the description of the high- 
pressure phase assemblages in the system CaAl^Si^Og, it is too complicated to be 
described by a single high pressure phase equation of state such as Jeanloz Sc 
Ahrens (1900). Although the shock temperature data can be roughly approxi- 
mated assunalng a single phase and choosing a temperature dependent model 
for Cy, the scatter of data around the best such single-phase curve in the 
pressure-temperature plane is still significantly greater than the estimated 



errors In the temperatare measurements. It Is more reasonable to interprert 
the shock-temperature data by assuming that a series of phase transitions con- 
tinue to occur along the Hugoniot. The data are consistent with three phase 
transitions— at about 55, 85, and 100 GPa— each with a transition ene/gy of about 
0.5 MJ/kg. The exact transition energies and pressures are not, however, v;ell 
constrained by the data. The most realistic model would include both a 
temperature-dependent Cy and phase transitions. It is also possible that when 
melting does occur, it does so incongruently, with a compositional difference 
between the melt and solid. Such a transition must be treated in the context of 
a multi-component system, in contrast to the stishovite melting transition, for 
which SiOg can be treated as a single-component system. A possible conse- 
quence of incongruent melting would be a spreading of the mixed phase region 
out over a larger pressure range, and a resulting shallow slope for the Hugoniot 
in the pressure-temperature plane. 

Success of modelling the time-dependence of the measured spectral radi- 
ances demonstrates the potential usefulness of the optical pyrometry method in 
studying optical properties of shocked material, in addition to measuring shock 
temperatures. By taking advantage of this fact the decay times of emitted light 
can be used to determine the absorption coefficient of the material behind the 
shock front, as in the case of anorthite glass, for which the absorption 
coefficient increases from 1 to >14 mm"^ as a function of shock-pressure (Table 
2-2). Future experiments can be designed to determine absorption coefficient 



as a function of shock-pressure, and the shock state can be further character- 
ized in this way. 
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Chapter 3 

Hugoniot Equation of State of Anorthite Glass and 
Lunar Anorthosite 

ABSTRACT 

Twenty-one Hugoniot experiments were conducted on an 
amorphous material of anorthite composition, in the pressure 
range from 8 to 120 GPa, using both routine and new methods. Two 
Hugoniot measurements at about 120 GPa were made on Lunar gab- 
broic anorthosite. Theoretical Hugoniots are constructed for both 
materials assuming they disproportionate to their component 
oxides. These accurately predict the P-p behavior of the lunar 
anorthosite Hugoniot at 120 GPa and the anorthite glass Hugoniot 
above 50 GPa, but overestimate the shock temperatures of anorth- 
ite glass. They do not predict the release paths of either material. 
We conclude that the mixed oxide model is a good description of 
the bulk properties of the high pressure phases of anorthite, but 
does not represent the actual phases. A significant enrichment of 
the earth’s lower mantle is not precluded by the bulk properties of 
the anorthite high pressure phases. 




-58- 


Introductioa 

Anorthite is a mineral of considerable geophysical importance because it is 
among the first phases to condense in nebula models of the formation of the 
solar system (Grossman ic Jjarimer, 1974). If the accretion process began before 
condensation was complete, the lower mantle would be expected to be enriched 
in refractories (Tureklan & Clark, 1969), Such an inhomogeneous model would 
predict that a high pressure phase assemblage derived from anorthite is a major 
component of the lower mantle. Calculations based on cosmochemlcaL abun- 
dances of Si, Mg, Al, and Ca (Ross & Aller, 1976) predict a mean mass fraction of 
0.12 and 0.16 for an anorthite-composition phase in a predominantly enstatite- 
or forsterite-composition mantle, respectively. 

Due to its geophyv^ical importance, anorthite has been the subject of several 
shock wave studies. McQueen et al. (1967) conducted Hugoniot experiments on 
Tahawus anorthosite at pressures from 15 to 95 GPa. The samples used in this 
study were crystals, which were 90% plagioclase with .49 mole fraction anorthite, 
and initial densities ranging from 2.70 to 2.79 Mg/m® . McQueen et al. (1967) 
inferred from their data that above about 33 GPa the Hugoniot enters a regime 
in which the shock states correspond to a high pressure phase with a zero pres- 
sure density of between 3.48 and 3.53 Mg/m®, and an initial isentropic bulk 
modulus of 88 to 112 GPa. This high pressure phase gives a close match to the 
seismologLcally-determined lower mantle density, but the inferred bulk modulus 
is considerably less than the lower mantle at comparable conditions. 


Additional Hugoniot experiments of Ahrens et al. (1969a) on plagioclase 
from 18 to 66 GPa also indicate a transition of this materia;' to a high-pressure 
high-density polymorph above 30 GPa. It was concluded in this study, and by 
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Ahrens et al. (1969b) that the high-pressure polymorph corresponds to the hol- 
landite structure. 

Jeanloz Sc Ahrens (1980) carried out a series of Hugoniot experiments on 
lunar anorthosite 60025, irhich is about 10% porous, and single-crystal Miyafce- 
zima anorthite. Shock pressures attained ranged from 40 to 120 GPa. The 
porous data were used to determine the Gruneisen parameter, and a thermal 
equation of state was constructed for the high pressure phase. The zero pres- 
sure density and initial isentropic bulk modulus were concluded to be 3.40 
Mg/m^ and 86.5 GPa, respectively for the calculated Hugoniot which gave the 
best fit to the data inferred by the •authors a§ belonging in the high pressure 
regime, above 33 GPa. This is in close agreement with the lower values of 
McQueen et al. (1967). Jeanloz Sc Ahrens (1980) also concluded that, despite a 
close match in density to the lower mantle, the high pressure phase or phase 
assemblage of anorthite is too compressible to constitute a major lower mantle 
component. 

Boslough et al. (1983a) conducted shock temperature experiments on 
anorthite glass in order to further study its high pressure thermal behavior at 
pressures from 50 to 115 GPa and temperatures of 2500 to 5600"K. They found 
significant deviations in measured temperatures from what would be calculated 
assuming a simple single-high-pressure-phase model in this pressure range, and 
concluded that as many as three phase transitions might occur between 50 emd 
115 GPa, at about 55, 85, and 100 GPa. Shock temperature experiments on 
anorthite glass were also carried out to about 40 GPa by Schmitt Sc Ahrens 
(1983). The high temperatures (m 3000®K) measured in these experiments were 
assumed to be at the melting point of the high pressure phase existing along 


shear bands in the shocked materiaL 


The purpose of this paper is to present the first Hugoniot data on anorthite 
glass, and discuss it in context with the previous porous and single crystal 
Hugoniot data and shock temperature data to generate a more complete self- 
consistent high pressure equation of state. For comparison, new Hugoniot data 
are presented for a lunar recrystallized and brecciated gabbroic anorthosite 
(Apollo 15,418). This rock was the subject of previous shock wave work by 
Ahrens el ai. (1973) at pressures up to 28 GPa. 

Experiment 

Several different experimental techniques were employed to determine the 
Hugoniot of anorthite glass, all based on measuring the shock velocity through 
the sample. Determination of partial and fully released states in some samples 
was made by the use of low density buffer material and in two cases the free sur- 
face velocity was measured using an inclined mirror affixed to the sample 
(Ahrens Sc Gregson, 1964). Experiments were performed using either a 40 mm- 
bore propellant gun (Ahrens et al., 1971) or a two-stage light-gas gun. 

First, the standard Hugoniot experiment, as described by Jeanloz Sc Ahrens 
(1977) was used (Fig. 3-1). Second, anorthite ’’rider" mirrors were used in lieu of 
fused quartz mirrors on Hugoniot experiments which were primarily for other 
samples (Fig. 3-1). Third, shock transit times were measured from the oscillo- 
grams used primarily to determine the shock temperatures (Fig. 3-2) and 
finally, transit times were determined by obtaining streak camera records of 
light transmitted from the sample via optic fibers in several of the shock tem- 
perature experiments (Figs. 3-3 and 3-4). 
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Fig. 3-1. Configuration of standard Hugoniot experiments including schematic 
representation of streak. Times to, ti, tg, and ts correspond to shock wave 
arrival at: driver-arrival mirror interface, arrival mirror free surface, sample- 
buffer mirror interface, buffer rnirror-bufler arrival mirror interface, respec- 
tively. In ideal experiment (no shock tilt or bowing), transit times are: (tg— tq) 
for sample, (ti-tq) for arrival mirror, and (tg-tg) for buffer mirror. In the 
"rider" experiment, arrival mirrors are made of anorthite glass. 
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Fig. 3-2. Oscilloscope record from 750 nm channel of optical pyrometer: shot 
LGG 132 (An6T) to 40 GPa. At time t^, shock wave arrives at driver-sample inter- 
face. At t|, shock reachej sample free surface. 


-81 - 


ORIGINAL* PA'.:''" tP 

OF POOR QUALfiT 



Fig. 3-3. Configuration of shock-temperature experiments in which Hugoniot 
data are also obtained. In one method, the shock transit time is taken from ths 
optical pyrometer oscilloscope records of the light radiated from the shocked 
sample. In the other method, light is collected by optical fibers and conducted 
to a streak camera. 
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Fig. 3-4. Streak camera record from shot LGG 136 to 92 GPa. illustrating use of 
optical fibers in conjunction with shock temperature experiment to determine 
shock transit time. At time to, shock wave enters sample from opaque driver, 
and fiber collects thermally radiated light. At t|, shock reaches sample-fiber 
interface and light intensity decreases. Time calibration is by a laser modulated 
at 60 V'Hz by a Pockels cell. Three of six fiber traces are shown. 


In the standard Hugoniot experiment (Fig. 3-1), a shock wave is generated 
by the impact of a light-gas gun or 40mm gun-launched flyer plate with a sta- 
tionary driver plate, on which are mounted Liie polished sample and two silvered 
quartz or polycarbonate (Lexan) arrival mirrors. After traveling through the 
driver plate, the shock enters the sample and mirrors. When the shock wave 
arrives at the driver-mirror interface, the lapped surface of the driver is driven 
into the specularly reflecting surface of the mirror, resulting in a rapid change 
in reflectivity which can be observed, by monitoring the reflected light with a 
high speed Model 339B Beckman-Whitely continuous writing rotating mirror or 
TRW model 21A image converter streak camera. The shock wave continues 
through the sample to its other surface, where it impinges on a fused quartz or 
Lexan buffer mirror which changes its reflectivity in the same manner, providing 
a measure of the shock transit time through the sample. When the shock wave 
reflects from the buffer-sample interface, a rarefaction to a lower pressure 
state propagates back into the anorthite. Simultaneously a shock wave with the 
same pressure and particle velocity is driven into the buffer material. The 
arrival of this lower amplitude shock wave at the buffer material free surface is 
again detected by the loss of reLectivity of a buffer arrival mirror, which pro- 
vides a measure of the shock transit time in the buffer. In some cases the buffer 
material was replaced by an inclined mirror whose angle of extinction allowed 
determination of the free surface velocity of the shocked sample. Solution of 
the Reimann integral (Rice et aiL., 1958, Lyzenga & Ahrens, 1978) allowed deter- 
mination of a partial (buffer) or full (inclined mirror) release state in the sam- 


The "rider*' method is a modification of this experiment, in which the sam- 
ple is not anorthite glass, but some other material for which the experiment is 
primarily intended (Fig, 3-1). The arrival mirrors, however, are fabrioated out of 
emorthite glass. When the shockwave arrives at the driver-sample interface the 
silvered surfaces of the anorthite arrival mirrors lose their reflectivity in the 
usual way. The top, unsilvered surface of the mirror material continues to 
reflect at about 5% this intensity, and does so until the shock wave reaches it. 
The light signal disappears within « 2 ns of shock arrival, either as a result of 
shock tilt misaligning the free surface of the mirror in such a way that the light 
reflected from the xenon light source misses the aperture of the streak camera, 
or due to disruption of the specular reflecting quality of the surface by an unk- 
nown process. This appears on the streak as a lighter image, from which the 
transit time through the anorthite glass mirror can be determined. The anorth- 
ite glass mirrors still serve the same functional purpose as the normal mirror in 
determining the Hugoniot of the primary material. The states determined for 
the anorthite glass by this method are not as precise as those determined by 
the standard method, because the rider mirrors are about half the thickness (« 
1.1 mm) of the standard samples. This doubles the fractional error in the shock 
velocity measurement. A release state is not determined for anorthite glass 
with this method. 

Hugoniot states were also determined in the shock temperature experi- 
ments. The most straightforward method was to simply take the shook transit 
time from the oscilloscope shot records (Fig. 3-2), which display light intensity 
as a function of time (Boslough ei al„ 1983a). This technique is less precise than 
the standard method for two reasons. First, the time resolution is limited by the 


rise times of the photodiodes, which are 1 to 7 ns. Second, the field of view of 
the photodiodes (2.3 to 8.4 mm) is much larger than the £ vs dth of the streak 
camera (0.13 mm), so if lae shock wave is not perfectly parallel to the sample 
surfaces (because of projectile tilt or slight misalignment) this will not be 
resolved. For sortie samples, the optical quality of the sample may preclude the 
use of this method in determining shock transit times (Boslough ei al. 1983b). 

Several of the shock-temperature targets were designed to allow the streak 
camera to be used in measuring the shock-transit time, with the intent on 
improving the precision over reading the oscillograms. Several 0.25 mm diame- 
ter, Poly-Optical Products Co„ acrylic optical fibers with polished ends were 
attached to the free surface of the samples, and held in place by the edge mask 
(Fig. 3-3). When th=! shock wave enters the sample, thermal radiation is col- 
lected by the optic fibers and carried to the streak camera. When the shock 
reaches the sample-fiber interface, the light intensity changes due to the 
impedance mismatch as the shock is driven into the fiber itself. The transit 
time is obtained by measuring the length of the bright streak (Fig. 3-4). This 
technique is not degraded by shock wave tilt, as the field of view of the fiber is 
on the order of the fiber diameter itself. The time resolution is decreased some- 
what, however, by the fact that the fiber diameter is 0.25 mm, as opposed to the 
slit width in the standard experiments of 0.13 mm. By using 0.13 mm diameter 
fibers, or mounting larger fibers behind the 0 |3 mm wide slit, this method can 
achieve the same precision as the d Hugoniot experiment. A shock tem- 

perature of about 3500“K in the g’ampl#, tjprregponding to a spectral radiance in 
the visible of about 10^^ Wsr~^ra''^^, Is rseOJired ta produce enough light intensity 
to give a sufficient exposivJt of tht I'tlfii ||1 ihe streak camera when 0.25 mm 
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diameter fibers are used. If finer fibers are used, the required shock tempera- 
ture be higher. The required spectral radiance is approximately inversely 
prupurliunul to the fiber diameter. 

Results 

The results of 21 Hugoniot measurements on anorthite glass are presented 
in Table 3-1 and shomm in Figs. 3-5 through 3-8. Previous data for Izu Island 
anorthite (Jeanloz &; Ahrens, 1980), porous lunar anorthosite 60025 (Jeanloz & 
Ahrens, 1978) and Tahawus anorthosite (McQueen et aL, 1987) are also shown in 
Fig. 3-6 for comparison. 

All but six data from the present study are in the high pressure phase 
region, and are plotted in the U,-Up plane in Fig. 3-7, high pressure phase 
data can be fitted by the linear relationship: 

U, = Co + sup (3-1) 

A least squares fit gives the values Co=1.64±0.26 km/ s and s=1.82±0.07. This 
equation can be used for impedance match solutions, and was used to determine 
the pressures in the shock-temperature experiments by Boslough et aL. (1983a). 

The remaining data lie in the mixed phase region or in the low pressure 
phase region and form an approximately horizontal trend in Ug-Up space (Fig. 
3-7). Two samples have anomalously high shock velocities relative to the 
remainder and were excluded from the diagram. Poor shot records made these 
two experiments difficult to interpret and although they are included in Table 3- 
1, they will not be discussed further. A similar horizontal trend in U,-Up space 
has been described by Wackerle (1982) for the mixed phase region of fused 
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Table 3-1 


Anorlhite Qaas HugonioL and Release Data 


Shot 

Flyer Driver Projectile 

Uaterial Material Velocity 
(km/s) 

Initial 

Density 

(Hg/ VD?) 

Shock 

Velocity 

(km/s) 

Particle 

Velocity 

(km/a) 

Pressure 

(CPa) 

Density 
(Hg/ m®) 

Method 

40-577 

Ai-3024 Al-8024 0.95 

± .01 

2.691 

.002 

6.85 
± ,12 

0.45 
± .01 

6.32 
± ,12 

2.681 

0.006 

(e) 


Partial Release State 1 
(Lexan huffer) 


3.48 
± .05 

0.73 
± .03 

3.05 

±.1B 

2.76 
± .03 

(a) 


Partial Release State 2 
(Graphite Foam buffer) 


1.60 

±.02 

0.62 
± .02 

1.06 

±.04 

2.65 
± .01 

(a) 

40-568 

Lexan Al-2024 1.97 

i: .01 

2.697 

J;.002 

6.42 
± .05 

0.464 
± .004 

8.37 

±.06 

2.912 
± .003 

(a) 


Partial Release State 
(Lexan buffer) 


3.46 

±0.03 

0.72 
± .02 

2.96 

±.11 

2.827 

±.014 

(a) 


Final Release State (e) 
(Enclined mirror) 


- 

1.02 
± .03 

0.00 

2.64 

±.04 

(a) 

40-581 

Al-2024 M-2024 1.00 

± .02 

2.693 
± .002 

6.36 
± .04 

0.49 
± .01 

8.46 

±.18 

2.920 
± .006 

(a) 


Partial E2elease State 
(Lexan buffer) 


3.50 
± .03 

0.745 
± ,019 

3.12 
± .11 

2.822 
± .021 

(a) 


Final Release State (e) 
(Laclined mirror) 


- 

0.87 
± .02 

0.00 

2.78 
± .04 

(a) 

40,580 

Al-2024 A1-2C«4 1.105 

± .010 

2.692 
± .002 

6.745 
± .005 

0.534 
± .006 

9.69 
± .10 

2.923 
± .003 

(a) 


Partial Release State 1 
(Lexan buffer) 


3.47 

±.04 

.726 
± .025 

3.01 
± .14 

2.876 
± .013 

(a) 


Partial Release State 2 
(Polystyrene foam buffer) 


1.07 
± .05 

0.74 
± .04 

0.045 
± .005 

2.686 

±.016 

(a) 


I 
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Table 3-1 (cont’d) 


Shot 

Flyer 

Driver 

Projectile 

Initial 

Shock 

Particle 

Pressure 

Density 

Method 


Uaterial 

Material 

Velocity 

Density 

Velocity Velocity 







(hm/s) 

(Hg/m®) 

(km/s) 

(km/s) 

(GPa) 

(Mg/ m°) 


40*567 

Al-2024 

Al-2024 

2.266 

2.693 

6.46 

1,179 

20.52 

3,295 

(«) 




:1:.003 

± .002 

±.06 

±.006 

± .11 

± .010 


Partial Release State 1 


3.53 

0.784 

3.227 

3.19 

(a) 


(Lexan biifler) 


±.04 

±.025 

± .144 

± .02 


Partial Release State 2 


1.04 

1.52 

0.165 

3.235 

(a) 


(Polystyrene foam buffer) 


±.02 

± .02 

± .004 

± .013 

40*566 

W 

Cu 

2.455 

2.693 

6.40 

2.279 

39.28 

4.162 

(a) 




0.005 

± .002 

± .04 

± .006 

± .21 

± .018 


Partial Release State 1 


6.24 

2.49 

16.59 

4.14 

(la) 


(Lexan buffer) 


± .08 

±.05 

± ,62 

± .03 


Partial Release State 2 


4.09 

0.15 

16.67 

3.66 

(a) 


(Graphite foam, buffer) 


± .05 

± .04 

± .37 

± .04 

LGG 132 

Cu 

Cu 

3.72 

2.690B 

6.508 

2.731 

47.8 

4.636 

(c) 




±.10 

± .0010 

±.062 

±.078 

± 1.5 

± .112 

LOG 131 

Cu 

Cu 

4.038 

2.6909 

7.027 

2.929 

55.4 

4.614 

(c) 

t 



± .002 

± .0013 

± .081 

±.008 

± 0.5 

± .047 

LLL 091660 

Cu 

Cu 

4.536 

2.6900 

7.717 

3.247 

67.4 

4.644 

(c) 




±.006 

± .0010 

± .103 

±.011 

± 0.7 

± .056 

LGG 136 

Cu 

Cu 

5.255 

2.6953 

8.187 

3.753 

82.B 

4,977 

fd) 




± .003 

± .0012 

± .240 

±.025 

± 1.9 

± .152 

LLL 091680 

Cu 

Cu 

5.455 

2.6896 

8.593 

3.867 

89.4 

4.BBB 

(o) 




± .010 

± .0011 

± .050 

±.009 

± 0.4 

± .030 

LGG 136 

Cu 

Cu 

5.506 

2.6923 

a769 

3.887 

91.8 

4.636 

(d) 




± .003 

± .0011 

± .085 

±.009 

± 0.7 

± .046 

LGG 137 

Cu 

Cu 

5.563 

2.6927 

6.782 

3.930 

92.9 

4.874 





±.006 

±.0011 

±.093 

±.011 

±0.8 

±.052 
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Shot 

Flyer 

Driver 

Projectile 

Initial 

Shock 

Particle 

Pressure 

Density 

Method 


Material 

Material 

Velocity 

Density 

Velocity Velocity 







(km/s) 

(Kg/m'») 

(km/s) 

(km/s) 

(GPa) 

(Hg/ m®) 


LGG 124 

Cu 

Cu 

5.709 

2.6894 

8.956 

4,025 

96,9 

4.885 

(a) 




dt.003 

±.0018 

±.068 

±.008 

±0.6 

±.037 


Partial Release State 


8.512 

4.461 

B3.7 

4.566 

(a) 


(Fused quartz bufier) 


±.095 

±.051 

±1.9 

±.116 

LOG 114 

Cu 

Cu 

5.692 

2.6812 

9.352 

4.124 

103,8 

4.614 

(b) 




f.005 

±.0018 

±.152 

±.016 

±1.3 

±.076 

LGG 133 

Cu 

Cu 

6.166 

2.6912 

9,535 

4,315 

110.7 

4.915 

(b) 




:i;.004 

±.0018 

±.090 

±.010 

±0.6 

±.047 

LLL 062360 

Ta 

Ta 

5.678 

2.6893 

9.525 

4,345 

111.3 

4.945 

(a) 




±mo 

±.0011 

±.100 

±.012 

±1.0 

±.052 

LGG 134 

Cu 

Cu 

6.328 

2,6812 

9.785 

4.412 

116.2 

4.801 

(b) 




±.004 

±.0008 

±.160 

±.017 

±1.5 

±.061 

LGG 139 

Ta 

Cu 

5.562 

2.6944 

9.839 

4.444 

117.8 

4,914 

(d) 




±.003 

±.0012 

±.089 

±.010 

±0.8 

±.045 

LGG 129 

Cu 

Cu 

8.547 

2.6924 

9.784 

4,580 

120.6 

5.062 

(a) 




±.004 

±.0020 

±.106 

±.012 

±1.0 

±.060 


Partial Release State 


9.328 

4.899 

100.7 

4.935 

(a) 


(Ftised quartz buffer) 


±.153 

±.062 

±3.3 

±,104 

LGG 143 

Ta 

Ta 

8.1 to 

2.6903 

9,043 

i.699 

124,4 

5,148 

(b) 




±.005 

±.0014 

±.136 

±.013 

±1.4 

±.077 


(a) Standard equation of state experiment. 

(b) Mirror "riders" - (density from mean of other shots). 

(c) Transit time from shock temperature oscUliSeope records. 

(d) Transit time from optic fibers to streak camera. 

(e) Free surface velocity listed as particle velocity. 



Pressure (GPa) 




Density (Mg/m^) 


Fig. 3-5. Anorthite glass Hugoniot data in the tiigh-pressure regime from Table 
3-1 plotted in the pressure-density plane, with theoretical mixed oxide Hugoniot 
based on parameters in Table 3-2. Data were obtained by four experimental 
methods: (S) Standard Hugoniot, (R) "Rider" mirror, (O) Oscilloscope, (F) Optical 
ftber. These methods are described in text. 
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Fig. 3-6. Anorthite glass Hugoniot data from this study, anorthite crystal and 
porous anorthosite data from Jeanloz & Ahrens (1900), and anorthite crystal 
data from McQueen et cti., (1987). Solid and dashed curves give single phase and 
mixed phase regions, respectively, based on shock temperature data (Boslough 
et al., 1983). Mixed oxide theoretical Hugoniot based on parameters in Table 3-2 
gives good agreement with data above 55 GPa. Grossularite-Kyanite-Stishovite 
(G-K-S) theoretical Hiigoniot does not agree with data anywhere. 
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Shock Velocity (km/s) 



Particle Velocity (km/s) 

Fig. 3-7. Hugoniot data for anorthlte glass shown in the shock velocity-particle 
velocity plane. Best straight-line fit to high-pressure phase regime data is 
shown: U,= 1.84+1.82up. Calculated zerorpretiisure longitudinal sound velocity 
(vp) and bulk sound velocity (vp) are indicated. 
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Fig. 3-9. Anorthite glass Hugoniot data with measured release states. Curves 
are Birch-Murnaghan isentropes constrained by Kj,,=60.4 GPa, K'o,=15-16 and 
K5,=:67.0, K'o,= 1 1.5-12.5. Release paths suggest that irrever.sible compression 
can occur in anorthite glass at low pressures and are consistent with the calcu- 
lated range of isentropes. 
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quartz, 

None of the data show any indication of a two wave structure. This suggests 
that the data lie below the HugonLot elastic limit (H.E.L) or sufficiently far above 
it that the elastic wave has been overridden. Particle velocities calculated from 
free surface velocities which were derived from inclined mirror experiments are 
covnapatible with particle velocities calculated by Impedance match. This is con- 
sistent with the observation of a single shock wave. 

For crystalline feldspars H.E.L. states are dependent upon driving shock 
pressure and vary upwards from 3.3 GPa for albite, 4.5 GPa for anorthosite 
(Ahrens k Gregson, 1964) and 3.5 GPa for microcline (Ahrens & Liu, 1973). For 
crystalline and fused quartz the lowest H.E.L states measured by Wackerle 
(1963) were 4.5 GPa. The states measured in this work for anorthite glass yield 
pressures on the order of 8 GPa. These pressures are considered to be too high 
to represent states below the H.E.L. and are suggested instead to be plastic 
deformational states well above the H.E.L. 

Release data from the low pressure phase regime and at the low density end 
of the mixed phase region show generally steep release paths (Figs. 3-B k 3-9). 
Buffer materials used were polystyrene foam (po - 0.056 lig/rn^, c^ = 0.343 
km/s, s = 1.118) and graphite foam (po = LOBl Mg/m^ Co = 0.79 km/s, s = 1.30). 
Inclined mirrors were also used to help constrain the release path in two experi- 
ments. At the high density end of the mixed phase region the release path sug- 
gests reversion to a very low density phase. Fluid-like behavior above the H.E.L. 
has been observed by Ahrens & Rosenberg (1968) for polycrystalline quartz and 
has been predicted theoretically by Zamyshlyaev ef al. (1983). 


Pressure (GPa) 




CaAl2Si20Q 

Hugoniot State 1 study 
4^ Portially Released > pi ■ 
State 


Hugoniot State 
^ Partially Released 
State 


Jeanioz and 
Ahrens (I960) 
Crystal 


Theofiitical Release Path 






Density (Mg/m'^) 


Fig. 3-9. Anorthite glass and crystal Hugoniot data, with measured release 
states. Theoretical mixed-oxide Hugoniot and corresponding theoretical release 
isentropes from experimental shock states are shown in the high pressure phase 
regime. The lowest and two highest states in the high pressure region demon- 
strate release paths which agree well with the theoretical paths. The two states 
between 90 and 100 GPa. however, release along significantly shallower paths. 
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ham preamire equation of state 

The low {pressure equation of state can be constrained using the two lowest 
pressure data. The H.E.L. has probably been exceeded for these experiments. 
At shock states higher than the H.E.L for anorthite glass it is expected that 
resistance to shear is lost and that it behaves as a fluid. This phenomenon has 
been observed in crystal and fused quartz (Wackerle, 1962), perthltic feldspar 
(Grady aL, 1975) and sapphire (Graham & Brooks, 1971). It has been sug- 
gested (Grady et ai., 1975) that the loss of shem' strength is related to partial 
naelting in shear bands behind the shock front. Shock induced shear instabili- 
ties have also been reported for anorthite glass by Schmitt & Ahrens (1983) pro- 
viding evidence that anorthite glass loses all resistance to shear in common with 
other ceramic materials. 

In the region above the H.E.L, with no shear resistance, the Hugoniot 
should differ little from the hydrostat. In addition, at low pressures, thermal 
effects will be small amd the Hugoniot will be well approximated by the principal 
isenlrope. 

No direct measurement of the bulk modulus of anorthite glass has been 
made but Young’s modulus of an An 7 o glass measured by Birch &c Bancroft (1942) 
is 89,9 GPa. A strong correlation between the product of Young’s modulus (E) 
and ionic packing density (Vj, and the observed bulk modulus (K) for a set of 30 
glasses of varying composition is observed by Makishima & MacKenzie (1975). 
This relation is given by K = 1.2VtE with calculated from ionic radii and the 
component atoms normalized to glass density (Makishima & MacKenzie, 1975). 
From this linear relation, the bulk modulus is estimated to be 60.4 GPa based on 
E=89.9 GPa, the value measured by Birch & Bancroft (1942). Poisson’s ratio for 
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these values of K and. E is 0.252. 

Further constraints may be placed on the bulk modulus of anorthite glass 
using recent ultrasonic data from a variety of silicate glasses. Meister et al. 
(1980) measured the elastic moduli of a range of silicate glass compositions up 
to 0.0 GPa and noted a good linear relation between bulk modulus and silica con- 
tent. For anorthite glass this relation would predict a bulk modulus of 67 GPa, 
and a Poisson’s ratio of 0.276, in reasonable agreement with the values calcu- 
lated using the relation of Makishima Sc MacKenzie. 

With these values of bulk modulus and an assumed Birch-Murnaghan equa- 
tion of state the pressure derivative of the bulk modulus (K’,) may be deter- 
mined from the Hugoniot data. It is assumed that at these pressures the princi- 
pal isentrope is coincid(§|jt with the Hugoniot, For Ko„=60.4 GPa, K'ogSai5--16, and 
for Kog=67 GPa, Keg'll 1.5—12.5 (Fig. 3-8). These values of K'o, suggest rapid 
stiffening of anorthite glass for moderate increases in pressure. In contrast, 
K',<3 for rock silicate glasses at pressures less than 0.8 GPa (Meister ef al., 
1980). Elastic parameters have been measured for fused quartz from atmos- 
pheric pressure to 13.3 GPa (Shroeder ei al,, 1981). Above about 3 GPa, K'g 
increases in a non-linear fashion and at 8 GPa, K'g»i7.5. Boslough and Ahrens 
(1983b) suggest that jadeite, in which Al is in octahedral coordination, is formed 
on the anorthosite Hugoniot in the range 6-10 GPa. Development of some 
octahedral aluminum on the anorthite glass Hugoniot could lead to substantial 
stiffening and the observed strong dependence of bulk modulus on pressure. 
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High pressure equatioa of state 

The shock temperature data of Boslough et al. (l9B3a) indicate that at Least 
three phase transitions occur above 50 GPa on the principal Hugoniot of anorth- 
ite glass. The temperature data were used to estimate where the phase boun- 
daries cross the Hugoniot, and this information was transformed to the P-p 
plane in Fig. 3-6, where single- and mixed-phase regions are delineated along the 
Hugoniot. These phase transitions are not evident from the P-p Hugoniot data, 
except perhaps between 80 and 90 GPa, where there appears to be a sudden 
increase in density with pressure. 

The highest pressure phases are only stable on the Hugoniot in a density 
range of about 0.1 Mg/m®, or about 5% of their density, which is the same order 
as the scatter in the data. This makes it impossible to carry out a meaningful 
finite strain fit to the data in the individual regions of stabilityt in which several 
of the parameters; zero pressure density, Kj,,, and K',,, are constrained indepen- 
dently. One alternative is to fix onf parameter and determine what values the 
others must have in order to fit the data. This method does not, however, lead 
to a unique solution. 

Another strategy which can be used is to choose a candidate phase or 
assemblage, a priori and determine its high pressure properties. It can then be 
compared to the data directly and accepted or rejected as a candidate on this 
basis. This is similar to the approach of Davies & Gaffney (1973). We used a 
mixed oxide model as the candidate high pressure assemblage, with CaO in the 
BS phase, AI 2 O 3 as corundum and SiOg as stishovLte. ' 

llie isentrope for the composite mixed oxide assemblage was constructed 
from the isentrope s of the individual oxides. Molar volumes as a function of 
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pressure were determined for each of the oxides using the Birch-Murnaghan 
equation and published values of K,,, and K’o, (Table 3-2). The values for B2 CaO 
were taken from floslougli ei al. (1983b), corundum from Anderson (1973) and 
stishovite from Lyzenga ef al. (1983). The molar volumes at a given pressure 
were added to give the molar volume of the composite, from which the 
pressure-density isentrope was obtained to 100 GPa. Ko,=227 GPa and Ko,'=5.1 
were determined for the composite by carrying out a Birch-Murnaghan least 
squares fit to the calculated isentrope (Vasiliou A: Ahrens, 1981). The zero pres- 
sure bulk modulus calculated in this way is comparable to the Voigt-Reuss-Hill 
value of 228 GPa, and lies within the Voigt and Reuss bounds of 212 and 244 GPa, 
respectively (Watt et al. 1976). 

The Gruneisen parameter ( 7 ) for the composite was calculated by taking 
the mean of the components, weighted according to mass fraction. Although 
this is not strictly correct (Duvall Sc Taylor, 1971) it is a reasonable approxima- 
tion. The energy of transformation between anorthite glass and the metastable 
high pressure assemblage at standard conditions (Et,) was 'btained from the 
difference between the sum of Etr for the oxides, and Etj. for anorthite glass, with 
the stable oxides as the reference state. 

A Hugoniot was calculated for the composite mixed oxides, with anorthite 
glass as the initial state. This theoretical Hugoniot is in agreement with shock 
wave data above 60 GPa (Figs. 3-5 A 3-6), but the temperatures predicted by this 
model are as much as 1500°K higher than those measured (Boslough et al., 
1983a), In order to fit data in both the P-p and P-T planes, a somewhat larger 
transition energy is required. 
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Table 3-2 

Mixed Oxide Parameters - AnorUute glas 


Po (g) 4.00 3,986 4.29 4,07 

Ko, (GPa) 121 250 306 222 

5,2 4,0 5,4 g 

^ 1-32 1.38 1,44 

n LQ 1.0 3.2 1.7 

Etr(kJ/g) 2.1 0,00 0.92W 0,87(*’) 

mi .202 .306 .432 1.00 

.216 .374 .410 1,00 

Reference (c) (d) (e) 

(a) Robie (1979) 

formation of anorthito gUsa oxides at standard 

(c) Boslough et al. (1983) 

(d) Anderson (1973) 

(e) Lyzenga et al. (1983) 
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For comparison, a Hugoniot was calculated in the same manner for an 
assemblage consisting of Ca 3 Al 2 (Si 04 )g (grossularite), AlgSiOg (kyanite), and SiOg 
(stishovite), which is the stable phase at lb GPa according to Liu (1978). This 
assemblage has a zero pressure density of 3.61 Mg/m^ and Voigt-Reuss-Hill bulk 
modulus of 170 GFa. A range of Hugoniots was calculated for this phase, with 3.0 

2.0 and 7=7o(^^* calculated Hugoniot is steeper than the observed 
P 

Hugoniot (Fig. 3-6) and crosses it at about 30 GPa. 

Discuanon 

It was concluded by Boslough et al. (1903a) that, due to the complicated 
behavior observed in the anorthite Hugoniot in the P-T plane, a simple equation 
of state is insufficient to describe this material at high pressure. It is evident, 
however, that the mixed oxide assemblage is a good model for the behavior in 
the P—p plane above 60 GPa. whereas the grossularite-kyanite-stishovite assem- 
blage is not a good description in any pressure regime. The Jeanloz- Ahrens 
(1900) equation of state adequately describes the data in both P-p and P-T 
planes from 33 to 60 GPa. 

Comparison of the calculated mixed oxide Hugoniot to data in the P-T plane 
is evidence that the high pressure phase is not mixed oxides, but a phase with 
similar bulk properties and a greater energy of transformation. Above 55 GPa, 
the CaAlgSigOa undergoes further phase transitions, and each higher pressure 
phase has a larger total energy of transformation. These higher phase transi- 
tions are sufficient to lower the shock temperature, but are not significant in the 
P-p plane, implying only a small volume change is associated with each phase 
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transition. The porous Hugoniot was determined experimentally by Jeanloz & 
Ahrens (1978, 1980), and does not agree with that calculated using the mixed 
oxide model, which predicts liigher densities at a given shock pressure. In light 
of the complexity of the high pressure properties of CaAlaSigOe, a reasonable 
explanation for this is that the porous data do not represent the same phase as 
the non-porous data, but rather a lower pressure or higher temperature, less 
dense phase. 

It must be remembered that, even though the mixed oxide Hugoniot fits the 
data well globally, the true Hugoniot goes through several phase transitions, and 
each phase is only stable along a short segment of the Hugoniot, where material 
properties such as the bulk modulus and Gruneisen parameter may differ 
greatly from those of the mixed oxides. Perhaps a better determination of the 
local compressibilities would be from partial release data. Fig. 3-9 shows that, in 
fact, the release paths do vary greatly with pressure along the glass and single 
?!:pystal Hugoniots. Of the four high pressure phase regimes inferred from the 
shock temperature data by Boslough et al. ( 1983a), three are sampled by the 
five highest pressure Hugoniot-release experiments in this study and that of 
Jeanloz 8c Ahrens (1980). Two experiments sample the highest pressure regime, 
above 100 GPa. The release paths are within experimental error of those 
predicted by the mixed-oxide model. The two experiments in which anorthite is 
shocked into the next highest phase regime, between about 89 and 96 GPa, show 
significantly different release behavior, indicating release paths much shallower 
than predicted. The region of stability inferred for this phase is so small and the 
Hugoniot data so scattered that is not possible to correlate these release paths 
to a measured Hugoniot. One experimental data point, at about 60 GPa, lies just 
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above the lowest inferred high pressure phase regime, and its release path is 
slightly steeper than predicted. Highly variable release paths were also 
observed in porous ariorthite (Jeanloz k Ahrens, 1978, 1980) and in some cases 
the release paths were steeper than the Hugoniot inferred in that study. Unfor- 
tunately, partial release data are generally of lower quality than Hugoniot data. 

It is interesting to compare the mixed oxide adiabat to that of the lower 
mantle, determined seismologically (Dziewonski, ef aX. 1975; Anderson & Hart, 
1976; Dziewonski k Anderson, 1981) in Figs. 3-10 through 3-12. In order to make 
a meaningful comparison, it is necessary to correct the theoretical mixed oxide 
adiabat to temperatures existing in the mantle. The geotherm of Stacey, (1977) 
was used. Whereas Jeanloz k Ahrens (1980) fomd the bulk modulus much too 
low for anorthite to be a significant component of the lower mantle, the mixed 
oxide node! gives a bulk modulus that is somewhat high, but closer to the 
seismologically determined lower mantle value. Thus, the high pressure phases 
of anorthite should not be precluded as major constituents of the lower mantle, 
especially in light of the fact that there are unresolved phases on the Hugoniot 
which have bulk moduli which might deviate significantly in either direction 
from that of the mixed oxide composite. 

Lunar anorthosite 

In addition to anorthite glass, two experiments were conducted on lunar 
gabbroic anorthosite (Apollo 15, 418). These experiments were carried out with 
the standard method (Fig. 3-1), and partial release states were determined in 
both cases. The results are presented in Table 3-3 and Fig. 3-13. This rock was 
modeled as .74 volume fraction anorthite (anga) and .26 volume fraction 
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Fig. 3-10. Adiabatic bulk modulus of the mixed oxide model as a function of 
pressure on principal adiabat and geotherm (Stacey, 1977). Also shown is bulk 
modulus of mantle as determined from seismological data (Dziewonski Se Ander- 
son. 1981), and anorthite adiabat of Jeanloz & Ahrens (1980). 
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Fig. 3-11. Bulk sound speed of the mixed oxide model as a function of pressure 
qn principal adiabat and geotherm (Stacey. 1977), compared to seismologicaily 
determined curve for the earths mantle (Dziewonski & Anderson, 1981). 
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Table 3-3 


Apollo 15418 HugonioL ami Release Data 


Shot 

Fljer DriTer Projectile 

Initial 

Shock 

Particle 

Pressure 

Density 


Ifaterial iiaterial Velocit}^ 

Density 

Velocity 

Velocity 




(km/s) 

iiHg/w?) 

(km/s) 

(Km/s) 

(GPa) 

(Mg/m“) 

LGG119 

T« T« 5,823 

2.8414 

9,811 

4.397 

122.6 

5,140 


±.005 

±.0014 

±,060 

± .007 

± 0.6 

± .031 


Partial Release State 


9.838 

4,752 

101,2 

4.996 


(IfaCl buffer) 


± .016 

± .014 

±0.6 

± .033 

LGG12S 

Ta Ta 5.830 

2.8014 

9.907 

4.380 

125.5 

5.183 


± .005 

±.0022 

±.035 

± ,005 

± 0.4 

± .016 


Partial Release State 


9.170 

8.446 

63.9 

3.613 


(Graphite foam buffer) 


±.031 

±,024 

±0.5 

± ,032 




Pressure (GPa) 
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Fig. 3-13. Lunar gabbroic anorthosite Hugoniot and release data. Mixed oxide 
theoretical Hugoniot based on parameters in Table 3-4 gives good agreement 
with data. One partial release state vras determined for each shock state, but 
because shock states are so close, the two partial release states outline an 
approximate release path. 
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enstatite (ena 4 ) by Ahrens et al, (1973), taken from a microprobe analysis. In 
order to construct a theoretical mixed oxide Hugoniot for this rock, as was done 
for unurLhile glass, we assumed the plagioclase fiaetiori to be pure aiiurthite, 
due to our ignorance of the high-pressure behavior of NagO, and the observation 
that anorthite and albite have similar Hugoniots (McQueen Bt al., 1976), The 
resulting mass fractions in the simplified model are ,70 for anorthite 
(CaAlgSigOa), .19 for enstatite (MgSiOg), and .11 for ferrosilite (FeSiOg). These are 
broken down into their five component oxides (Table 3-4), and a high pressure 
adiabat and Hugoniot were constructed in thjB sani© .manner as for anorthite 
glass, The theoretical Hugoniot constructs^d in this way is seen in Fig. 3-13 to 
give excellent agreement to the Hugoniot data, and to the initial slope of the 
measured release path. As m the case of pure anorthite, it is not likely that the 
lunar rock disproportionates into its component oxides at these pressures, but 
the mixed oxide model gives a good description of the bulk properties of what- 
ever phases exist there. The second release datum indicates a large density 
decrease upon unloading, consistent with shock-vaporization of material. 
Entropy-gain calculations of Ahrens & O'Keefe (1977) and of Boslough &; Ahrens 
(19B3a) indicate that incipient vaporization may occur at shock pressures above 
102 GPa for gabbroic anorthosite and above 92 GPa for anorthite, respectively. 
However, this is not observed in the release paths of shocked anorthite glass. 


Table 3-4 

Mixed Oxide Paramelers - Lutiaf Gabbroio Anorthosite 



References 


(a) Robie (1979) 


(b) includes .35kJ/g energy of formation of anorthite and .35kJ/g energy of for- 
mation of enstatite from oxides (Robie, 1979) 


(c) Boslough et al. (19B3) 

(d) Anderson (1976) 

(e) Ly/enga et al. (1973) 

(0 Va Esiliou and Ahrens (1981) 
(g) Jeanioz and Ahrens (19BU) 
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Concliuiona 

The experimental Hugoniots of anorfchite glass above 60 GPa and lunar gab- 
broic anorthosite at 120 GPa are well described by the Hugoniots calculated for 
their composite oxides, based on the individual properties of the oxides. The 
actual phases which are stable along the HugonLot of anorthite, however, have 
larger energies of transformation. The zero pressure density of the mixed-oxide 
high pressure phase assemblage model is 4.07 Mg/m^, and the zero-pressure 
bulk moduluu is 228 GPa. These values are comparable to those for the lower 
mantle, indicating that CaAlsSigOg, or its high pressure reaction products, can- 
not be precluded as major components of the lower mantle on the basis of its 
equation of state. Such a refractory component could therefore be present in 
the lower mantle with a mass fraction in the range 0.12 to 0.16 as predicted by 
cosmochemical abundance considerations (Ross Sc Aller, 1976), or in greater 
concentrations as expected by inhomogeneous accretion models (Turekian k 
Clark, 1969). 
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Chapter 4 

Shock Wave Properties of Anorthosite and Gabbro 

ABSTRACT 

Shock wave experiments have been conducted in San Gabriel 
anorthosite to peak stresses between 5 and 11 GPa using a 40 mm- 
bore propellant gun. Particle velocities were measured directly at 
several points in each target by means of electromagnetic gauges. 
Hugoniot states were calculated by determining shock-transit 
times from the gauge records. Measured sound velocities indicate 
a loss of shear strength upon shock compression for both rocks, 
with the strength loss persisting upon release to zero stress. 
Sound velocities of anorthosite shocked to peak stresses between 6 
and 10 GPa were measured to be between 5.1 and 5.3 upon release 
to zero stress. Stress-density release paths in the anorthosite indi- 
cate possible transformation of albite to jadeite + (quartz or 
coesite), with the amount of material transformed increasing from 
about .4 to .7 mass fraction as the shock stress increased from 6 to 
10 GPa. Electrical interference effects precluded the determina- 
tion of acurate release paths for San Marcos gabbro. Because of 
the loss of shear strength in the shocked state, the plastic 
behavior exibited by anorthosite indicates that calculations of 
energy partitioning due to impact onto planetary surfaces based 


-98- 


: 


on elastic-plastic models may underestimate the amount of inter- 
nal energy deposited in the impacted' surface material. 


Introduction 

The plagioclase feldspar-bearing rocks anorthosite and gabbro are impor- 
tant components of the lunar and terrestrial crusts. It is necessary to under- 
stand the behavior of such rocks under high dynamic stress in order to model 
cratering processes which result from hypervelocity impacts, ami to character- 
ize the stress history of rocks which have been subject to shock loading on 
planetary surfaces such as the moon and in some meteorites. Shock-wave stu- 
dies of these and similar materials have been conducted in the past (Ahrens et 
al., 1989, McQueen et aL, 1967, Jeanloz and Ahrens, 1980, fioslough et al., 1983a), 
but in these studies the data are limited to the Hugoniot state and in some cases 
a single state on the release isentrope. By employing particle-velocity gauges, a 
complete stress-strain history subsequent to shock-compression can be deter- 
mined, along with sound velocity information (Fowles and Williams, 1970, 
Cowperthwaite and Williams. 1971, Seaman, 1974). Particle velocity experiments 
supply detailed release paths, which provide better constraint for mechanical 
properties and polymorphism than is available with Hugoniot experiments edone. 

Peterson et al., (1970) used particle velocity gauges to determine release 
paths of playa alluvium, tonalite, and novaculite shocked to stresses up to 5 GPa. 
They attributed high rarefaction velocities and steep release paths in the 
stress-density plane to irreversible compaction. Grady et al., (1974) carried out 
experiments on polycrystalline quartz (novaculite) to 40 GPa using a 
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combination of particle velocity and manganin stress gauges to determine 
release adiabats and concluded that a partial quartz-stishovite transformation 
takes place above 15 GPa, with the quantity of material transformed an increas- 
ing function of peak stress, and that the Hugoniot states are not on the quartz- 
stishovite coexistence curve. Similar experiments were conducted on polycry- 
staUine quartz and perthitic feldspar by Grady et ai., (1975) and Grady and Murri 
(1976), who used manganin stress gauges to determine Hugoniot sound veloci- 
ties, and found that these rocks lose shear strength when shocked to pressures 
above 20 GPa. Larson and Anderson (1979) used particle velocity gauges to 
study limestone and tufif at lower stress levels (4 GPa), and attributed the 
observed time-dependent behavior to the closing of pores in these rocks. 

In this paper we present new Hugoniot data on San Gabriel anorthosite and 
San Marcos gabbro to 11 GPa. Release paths in the stress-density plane and 
sound velocities are reported, as determined from particle velocity data. 

Experimental Methods 

San Gabriel anorthosite samples were collected in the San Gabriel Moun- 
tains near Pasadena, California. This rock is highly variable in composition and 
texture and has been studied in detail by Carter (1982). The particular speci- 
men used in these experiments had randomly oriented plagioclase crystals with 
a mean grain size of 1-2 mm. Significant alteration was observed at grain boun- 
daries, and the composition is given in Table 4-1; determined with a petro- 
graphic microscope. An electron microprobe analysis indicated a composition 
of Angg for the plagioclase component. 
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Table 4-1 

San Gabriel AnorUiasLLe ComposiUon 
Mineral Volunae % 


Plagioolase * 

91 

White Mica 

6 

Epidate 

3 

Opaques 

trace 

Quartz 

trace 

Apatite 

trace 


kn^Q 

(Analysis by S. Rigdcn) 




The San Marcos gabbro was obtained near Escondido, California, and has 
been studied petrologic ally by Miller (1937). Samples from the same specimen 
used in Lhia study were used in impact and in spall-strength experiments by 
Lange et al„ (1983). The composition is tabulated in Table 4-3, Hock samples 
were cut into 1.5 mm thick, 3.2 by 4,5 cm rectangular slab's, and were bonded 
together with epoxy^ with a U-shaped copper particle-velocity gauge at each 
interface and one at the free surface. Gauges were photo-etched from 10 /im 
thick copper foil with a 12 fj.m polyamide (kapton) film backing. Targets were 
shock-loaded by impact of flat-faced polycarbonate (Lexan) projectiles fired 
from a 40-mm bore propellant gun at velocities from 1.4 to 2.4 km/s. Projectile 
velocities were determined using a laser intervalometer (Ahrens et al,> 1971). 

A uniform magnetic field at right angles to each gauge and to its direction 
of motion was supplied by a set of Helmholtz coils. An electromotive force V(t) 
was induced across the gauge element due to its motion through the field 

V(t) = BLUp(t) (4-1) 

where B is the magnetic field («1.9 kG), L is the effective gauge length (wl cm), 
and Up(t) is the gauge velocity: eq'!’ial to the particle velocity of the surrounding 
medium. These signals were recorded by an array of cathode ray oscilloscopes. 
The experiment is shown schematically in Fig. 4-1 and described in detail in 
Appendix IV. 

The geometry and time history of a typical experiment is illustrated by 
means of an x-t diagram in Fig. 4-2. The stationary rock target lies to the right 
of the origin, with four particle velocity gauges initially at intervals of 1.5 mm. 
The polycarbonate projectile approaches from the left and strikes the x=:0 


Mineral 


Plagioclase 
Amphlbole 
Cliriopyroxene 
Orthopyroxene 
Quartz 
Blotlte 
Opaques 
Alakali feldspar 
Calcite 
Chlorite 
Apatite 








Fig. 4“1. Schematic* drawing of particle velocity experiment, with major com- 
ponents indicated; A. Polycarbonate projectile, B. 40 mm-bore gun barrel, C. 
Timing laser, D. Photodectector, E. High-power switch (ignition), F. Capacitor 
bank. G. Helmholtz coils, H. Rock target, I. Self-shorting trigger pins, J. Fiducial 
pulse generator. K. Copper foil particle velocity gauge elements. 
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Fig. 4-2. Particle velocity experiment represented by x-t diagram. Projectile 
approaches stationary target from left and impacts at t::0. 
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surface of the target at time t=0, driving a shock wave to the right into the rock 
and to the left into the projectile. Each gauge is stationary until overtaken from 
tiiB left by the shuck wave, at which Lime it begins moving mth the particle velo- 
city associated with the Hugoniot state. The shock wave reflects from the free 
surface as a rarefaction wave, £uid each gauge again accelerates to the right as 
this wave passes through it from the right. 

Results 

Typical particle velocity records for shocked anorthosite and gabbro are 
shown in Fig. 4-3, The sudden increase corresponds to shock wave arrival, and 
the second increase corresponds to the acceleration from free-surface rarefac- 
tion. The gabbro records were found to be significantly noisier than the anortho- 
site records in all cases. This was presumed to result &om the presence of 
piezoelectric queurtz grains in the gabbro, whereas the anorthosite was relatively 
free of quartz. 

Digitized oscilloscope records for anorthosite shocked to 10 GPa are shown 
in Fig. 4-4. Synchronization of the four signals in time was achieved by means of 
a fiducial pulse received simultaneously at each oscilloscope, Shock-transit 
times were taken from the interval between arrival at different gauges and were 
used to determine shock velocity. The known projectile velocity and polycar- 
bonate Hugoniot (Marsh, 19B0) were used with an impedance-matph solution 
(McQueen et al., 1970) to determine the Hugoniot states achieved in the rocks. 
Eulerian sound velocities were determined from the transit time of the free- 
surface rarefaction front and the Hugoniot density. Hugoniot states and sound 
velocities for both rocks are given in Tables 4-3 and 4-4. 
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Fig. 4-3. Oscillograms of particle velocity-time profiles in (a) anorthosite and (b) 
gabbro. Gabbro signals are significantly noisier, presumably due to presence of 
quartz grains. 
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Fig. 4-4. Digitized particle velocity gauge records from experiment 40-571 in 
which anorthite was shocked to 10 GPa. 
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Table 4-3 

San Gabriel Anorthosite Shouk Wave Data 



Projectile 

Initial 

Shock 

Particle 

Pressure 

Density 

Sound 

Shot 

Velocity 

(km/s) 

Density 

(Mg/ma) 

Velocity 

(km/s) 

Velocity 

(km/s) 

(GPa) 

(Mg/ma) 

Velocit’ 

(km/s) 

40-572 

1.575 

2.656 

5.648 

.394 

5.91 

2.855 

6.90 


1.050 

1.013 

1.055 

1.022 

1.31 

1.011 

1.36 

40-570 

1.881 

2.653 

5.725 

.493 

7.47 

2.903 

7.59 


1.050 

X.007 

1.117 

1.020 

1.28 

1.014 

1.10 

40-571 

2.351 

2.653 

5.985 

.644 

10.23 

2.973 

8.39 


1.050 

1.007 

1.032 

1.020 

1.30 

1.011 

1.44 
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Table 4-4 


San Marcos Gabbro Shock Wave Data 


Shot 

Projectile 

Velocity 

(km/s) 

Initial 

Density 

(Mg/m^) 

Shock 

Velocity 

(km/s) 

Particle 

Velocity 

(km/s) 

Pressure 

(GPa) 

Density 

(Mg/m'^) 

Sound 
Velocity 
(km/ s) 

40-569 

1.594 

2.907 

6.487 

.286 

5.39 

3.041 

6.5 


i:.009 

±.008 

±.035 

±.011 

±.20 

±.004 

±.l 

40-573 

1.968 

2.886 

6.11 

.474 

8.36 

3.129 

8.5 


±.008 

±.001 

±.30 

±.016 

±.14 

±.022 

±.4 

40-555 

2.187 

2.892 

6.66 

.515 

9.91 

3.135 

— 


±.050 

±.007 

±.19 

±.020 

±.34 

±.015 


40-556 

2.242 

2.929 

6.898 

.515 

10.41 

3.165 

7.0 


±.032 

±.039 

±.042 

±.048 

±.0.84 

±.020 

±.5 

40-574 

2.418 

2.874 

6.404 

.608 

11.19 

3.176 

6.7 


±.025 

±.049 

±.016 

±.052 

±.77 

±.026 

±.5 
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The observed release waves are nonsteady simple waves and can be Inverted 
to stress-density release paths by numerically integrating the equations for con- 
servation of mass and linear momentum (Cowperthwaite and Williams, 1971) 


PoC(Up) 


(4-2) 


da 

0Up 


h 


Pot^(Up) 


(4-3) 


where p is the density, is the initial density, a is the stress, Up is the particle 
velocity, and h is the Lagrangian space coordinate along the direction of wave 
propagation. The Lagrangian sound velocity is determined by the finite 
difference approximation 


C(up) w 


At 


(4-4) 


where Ah is the initial distance between gauges and At is the transit time for a 
disturbance with particle velocity Up. In the case of anorthosite, the release 
paths are smooth and single-valued (Fig. 4-4). The oscilloscope records were 
digitized and integrated directly. The resulting release paths are consistent, 
and are plotted in Fig. 4-5. Because of the noisiness of the gabbro records (Fig. 
4-3b), it was necessary to approximate the rarefaction waves by smooth curves 
before integrating. Complete release paths were not obtained (Fig. 4-8), and the 
partial release paths are less consistent than those of anorthosite. Eulerian 

sound speeds, equal to ^^(u_), were also calculated for the release paths in 

P 

both rocks. These are plotted as a function of stress in Figs. 4-7 and 4-8. 


Stress (GPa) 


San Gabriel Anorthosite 

Hi 

Release path . 

CaAl 2 SipO 0 Glass 

(Boslough et a!., 1983) 

p Hugoniot state 
p Partially released 


Density (Mg/m'^) 


Fig. 4-5. Hugoniot states and release paths of San Gabriel anorthosite. Included 
are two Hugoniot states of anorthite glass (Boslough et oL 19B3). with respec- 
tive partial release states. 
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Density (Mg/m^) 


Fig. 4-6. Hugoniot states and partial release paths of San Marcos gabbro. Noisy 
particle velocity records precluded determination of release to zero stress. 
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DiacuMion 

Sound velocities of anorthosite from this study can be compared «o those 
measured ultrasonically for anorthosite rocks of aimilar composition (Birch, 
I960, Birch, 1961, Simmons, 1984, Anderson and Liebermann, 1968, Liebermann 
and RingTVOod, 1976). At T = 25®C and P > 0.4 GPa, longitudinal velocities have 
been measured in anorthosite vrith composition between An 4 g and Angg in the 
range 8.78 a: Vp ;a 7.47 km/s, and shear velocities in the range 3.97 s: v,sc 4.09 
km/s. Bulk sound velocities, 



are therefore in the range 5.04 ;S v^ :S 5.99 km/s. Fig. 4-8 demonstrates that in 
all three shock wave experiments, the sound velocity approaches 5.2 km/s as 
the rock releases to zero pressure, indicating a loss of shear strength upon 
shock compression which is never regained. This is analogous to observations by 
Grady et al.. (1975) for polycrystalUne quartz, where the loss of strength in that 
material was attributed to heterogeneous partitioning of shock energy into 
"shear bands" in Vfhich melting occurs. Optical observations consistent with this 
theory have been made for various minerals by Kondo and Ahrens (1983), and 
for anorthite glass by Schmitt and Ahrens (1983). 

Not as clear a case can be made for the loss of shear strength in gabbro, 
due to the lower quality of the data. However, in most of the experiments the 
measured sound velocities are consistently lower for the gabbro than for the 


anorthosite. 
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Release paths for anorthosite, illustrated in Fig. 4>5, are significantly 
steeper than the Hugoniot. Similar release behavior was observed in anorthite 
glass {initial density = 2.69 Mg/ni^ ) by Boslough ei aL. (1983b). Tiie reason Cor 
the denslfication of the glass upon release was attributed to irreversible <ioni- 
paction of the amorphous material behind the shock wave to a denser for.n. 
Partial release paths can be extrapolated to the anorthite single crystal density 
of 2.76 Mg/m®, This argument cannot be used to explain the release behavior of 
anorthosite, hovrsv^r. The phase transformation; 

NaAlSisOg (albite) -» NaAlSigOa (jadeite) +• SlOg (q^uartz or coesite) is a more rea- 
sonable explanation, particularly in light of the identification of shock-induced 
formation of jadeite from oligoclase in material from the Ries Crater In Germany 
(James, 1969). The assemblage (anorthite) 3 n + (jadeite 4- quartz)e 2 has a density 
of 3.01 Mg/m^, and (anorthlte) 3 s + (jadeite ^ coesite)gg has a density of 3.06 
Mg/ m^. At 293“^ NaAlSigOa + SiOg are the stable phases above 0.8 GPa (Clark, 
1866). The present experiments are well above the transformation pressure, so 
the only barriers to the phase transition are kinetic. The zero pressure density 
of the released anorthosite can be used to estimate the amount of material 
transformed. For example, the final densities are iji the range 2.72 to 2.79 
Mg/ m®, which requires the transformation of from 42 to 73% of the albite to high 
pressure phases. The possibility of pore-space crushing and irreversable com- 
paction would bring this estimate down somewhat. It is clear from the data that 
the fraction of albite transformed to high pressure phases is an increasing func- 
tion of peak stress. This result is analogous to the conclusion of Grady et ai., 
(1974) for higher pressures that the release path of shocked poly crystalline 
quartz is controlled by the quartz -* stishovite transition. The release data for 


gabbro, plotted in Fig, 4-7, are not of good enough quality to resolve whether the 
release paths are above or below the Hugoniot. 

ConcliuioM 

San Gabriel anorthosite loses shear strength upon shock to pressures 
greater than 6 GPa. Because strength effects can probably be neglected in this 
material, its release beha'rior can be attributed to the phase transformation of 
albite to jadite and quartz or coeslte. The amount of material transformed is 
controlled by kinetics, and is an increasing function of !fhock pressure. This 
may be due to the localization of shock heating into shear bands, in which the 
temperature is high enough for the phase transition to occur. The mechanical 
result of the shear banding is plastic beha\dor upon release. 

In calculations by O'Keefe and Ahrens (1977) of energy partitioning of a 
hypervelocity impact onto the surface of a planet, the release behavior of rocks 
was assumed to be elastic-plastic (Fowles, 1960, Davison and Graham, 1979). 
This assumption over-estimates the attenuation of the shock wave due to catch- 
up of the rarefaction wave from the free surface of the impacting body, because 
the assumed sound velocities are hig, sr than they actually are for plastic 
behavior, as determined in this study. Thus a larger quantity of surface 
material is shocked above a given pressure, and estimates by O'Keefe and 
Ahrens (1977) of the fraction of meteoroid kinetic energy deposited in the plane- 
tary surface material will be too low. 


H,J» 
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Appendix 1 

Wide-Band Optical Pyrometer Design 


The wide-band optical pyrometer is the temperature-measuring portion 
of the shock temperature-equation of state combination experiment at Caltech. 
It was designed with the intention of giving up some of the spectral- and time- 
resolution available with the six-channel narrow-band pyrometer (Lyzenga and 
Ahrens, 1979) in exchange for increased sensitivity. It is over two orders of mag- 
nitude more sensitive in spectral radiance than the narrow-band pyrometer, and 
can be used to measure temperatures lower than 2500'’K without 
preamplification of the signals. At the same time, the loss of spectral- and 
temporal-resolution appears to be negligible in terms of data quality. Sensa- 
llvily might In fad he further Increased by using even wider-hand filters. 

The description of the optical pyrometer will proceed In the direction of 
light flow, beginning with the target assembly {Figs. 1-1, 3-3). Samples are trans- 
parent solids about 17 mm in lateral dimension and about 3 mm thick. Both 
surfaces are lapped to optical quality, and the samples are glued to flat metallic 
driver plates with epoxy. Great care is taken to minimize the gap between 
driver plate and target, and it is generally on the order of one micron according 
to measurements made by counting visible fringes with an. optical flat. In some 
cases the contact surface of the sample is silvered before attaching to the driver 
plate in order to mask any flash which might occur when the shock wave transits 
the gap. 
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Fig. 1-1. Photograph of shock-temperature target assembly from direction of 
pyrometer. Aperture (8.35 mm diameter) in edge mask can be seen in 
reflection in turning mirror. Optical fibers carry light pulse to streak camera. 
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The driver plates are the same as those used in standard equation of 
state experiments. They are circular disks 35 mm in diameter and 2.0 mm 
thick, and are mounted on the same target assembly, in the aatue manner as 
the standard experiment. As in the standard experiment, shorting pins are used 
to generate a triggering pulse. Two pins are mounted on the driver plate 9 mm 
above and below the center point so as not to block light radiated from the sam- 
ple. The pins protrude through holes drilled in the driver plate, so the projectile 
reaches them about 200 ns before it collides with the target. The pins act as 
high speed switches which short upon impact. They are biased at 335 V, and 
when they close, a pulse is sent to the triggering circuits of the oscilloscopes 
and streak camera. 

Because the field of view of the pyrometer is greater than the diameter of 
the saumple, it is necessary that the edges be masked. It is also important that 
the outer portions of the sample, where edge eflects (two-dimensional flow) 
might occur must also be masked from view. For this reason an aluminum 
edge-mask is attached to the free surface of the sample. A circular hole, 6.35 
mm in diameter, and concentric with the axis of the driver plate (which is 
aligned with the gun barrel) is the only path by which light can get to the pyrom- 
eter. The edge mask also serves to hold the optical fibers— which carry light to 
the streak camera and make up the equation of state portion of the combined 
experiment— in position. 

The light generated in the sample is reflected by an aluminized mirror, 
turning it 90” toward the optical pyrometer port. The same mirrors and mirror 
mounts as used in the standard experiment are employed, but light is reflected 
toward the pjmometer, which is in the opposite direction from the streak 
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camera. Also, the mirror mounts are slotted to provide an extra degree of free- 
dom tor optical alignment. 

The port through which the light is reflected is covered with a glass win- 
dow which separates the evacuated target chamber from the external instru- 
ments, and protects the objective lenses from flying debris generated by the 
impact. The objective is made up of two identical sets of achromats, with focal 
lengths of 49.5 cm and clear apertures of 10.6 cm. The target is near the focal 
point of the first achromat, with the beam folded by the turning rrurror. 

The converging light is separated into four legs by three beamsplitters 
(Figs. 1-2, 1-3). The primary beamsplitter is an Ealing 22-0924 pellicle, or 
stretched-membrane. The secondary beamsplitters are identical Melles-Griot 03 
BTF 007 silvered plate-glass beamsplitters. Each of the four resulting beams 
falls upon a detector assembly , which consists of an interference filter, a focus- 
ing lens, and a silicon photodiode with a bias circuit (Fig. 1-4). The filters are 
25.4 mm in diameter and have wavelengths in the visible or near infrared, 
between the effective bandwidth limits of the silicon photodiodes of 400 and 1100 
nm. Most of the filters used have half-height bandwidths of about 40 nm. 

Each channel utilizes an RCA C30822 N-type silicon pin photodiode, with 
an active area of .20 cm®. In order to ensure that the image of the sample (and 
calibration lamp, which is larger), falls entirely within the photodiode active 
area, it was necessary to employ a glass focusing lens, with diameter 18 mm and 
focal Length 15 mm. This lens is positioned between the filter and photodiode, 
and demagnifles the image by a factor of about 0.2. The detector assembly is 
mounted on a three-dimensional translation stage to allow alignment and focus. 
The focusing lens and filter holder can be moved independently to permit 
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Fig. 1-2. Schematic of shock temperatiire-equation of state experiment. 
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Fig. 1-3. Photograph of four-channel wide-batnd optical pyrometer. 
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separate adjustment of focus and demagnlflcation. 

This optical configuration ensures that each detector views precisely the 
same image, only at a different wavelength. Thus, if the i-udiatlon is heterogene- 
ously emitted, as would be the case If shear bands or hot spots developed, the 
average brightness would be measured by each detector. This eliminates the 
possibility ol one detector viewing a hot spot while another views between hot 
spots. Also, by viewing the entire unmasked portion of the target, the sensitivity 
of the instrument is maximized (Fig 1-5). Because the image is demagnifled 
there is sufficient leeway in alignment and focusing; as long as the entire image 
falls on the active area of the detector, it doesn’t have to be centered or pre- 
cisely focused. 

The photodiode circuit is shown schematically in Fig. 1-4. The circuit sup- 
plies 45 V reverse bias from a battery to the photodiode. The signal is carried to 
oscilloscopes by a 50 0 cosasial cable, where it is terminated by a 50 Q termina- 
tor. The principal recording oscilloscopes are Tektronix 485 models. These are 
triggered by a pulse generated by the closing of shorting puis mounted on the 
target as described above. Backup recording is provided by Tektronix 475, 454, 
and 545 osciiloscopes. These are set at lower a'erisilivities and slower sweep 
rates, and are triggered by the same pulse that triggers the second projectile- 
velocity-determining flash x-ray. Tektronix series C-30 cameras with ASA 20,000 
Polai'oid type 612 ultra high speed instrumeiit recording film are used to record 
the oscillograms. 

The light gas gun on which these experiments are carried out is 
described by Jeanioz and Ahrens (1977). Details on the operation of the light gas 
gun and determination of the projectile velocity are discussed there. 
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Fig. 1-5. Diagram, shovring approximate sensitivity range of optical pyrometer. 
Curves give spectral radiance as a function of wavelength for a 6.35 mm diame- 
tet blackbody emitter at the position of the target. Dashed lines represent 
wavelength of detectors in typical experiment. Unstippled area is within sensi- 
tivity range of silicon photodiodes. 
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Appcndix n 

Wideband Optical Pyrometer Calibration 

The calibration procedure of the wide-band optical pyrometer was designed 
to minimize calibration error. By calibrating the instrument in situ, geometry- 
induced errors are u kiimized. The same oscilloscopes used to record the sig- 
nals are used in the calibrations, minimizing the number of oscillograms to be 
iv. 4 d, which reduces possible reading and transfer error. Finally, because a 
separate calibration is performed before each experiment, the effect of variabil- 
ity in replacement windows and turning mirrors is taken into account. Also, any 
contamination of optical components, such as dust or residue, becomes less 
important, and components can be replaced readily. 

A standard of spectral irradiance is placed in the position of the target in 
the impact chamber (Fig. II-l), and the pyrometer is aligned. Spectral irradi- 
ance is the power radiated per unit wavelength per unit solid angle, regardless 
of the surface area of the radiator. The standard source is a G.E. type 
Q6.6A/T4Q/1CL 200 Watt tungsten-filament quartz-halogen lamp, calibrated by 
Optronics Laboratories, Inc. It is mounted on a dummy target assembly and 
aligned in the same manner as a real target. It is driven by 6.50 Amperes from a 
Hewlett-Packard Model 6268B DC power supply constant current power supply. 
Current is monitored by a Keithley Model 173 Digital Multimeter, which meas- 
ures the voltage across a ,005 fl resister in series with the lamp. The light is 
physically choppeu at about 200 Hz, and the resulting square wave is recorded. 
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I Fig. 11-1. Photograph of pyrometer calibration set-up. with standard lamp in tar 

I get position and chopper in place. 
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The amplitude of the wave is divided by the spectral irradiance of the lamp at 
the wavelength at which it is filtered, to determine the calibration factor. This 
calibration factor is multiplied by the amplitude of the shock record to get the 
spectral irradiance of the shocked material. To determine the spectral radi- 
ance, from which the temperatiire is calculated, the spectral irradiance is 
divided by the area of the hole in the edge mask. The solid angle needs not be 
known, as long as the calibration lamp and the target are in the same position. 

The measured spectral irradiance in the wavelength range of interest of the 
calibration lamp is shown in Table II- 1, as supplied by Optronics Laboratories. 
Note that the units are /i,W/cm®nm at 50 cm, where the cm® at 50 cm is a meas- 
ure of solid angle and can be converted to steradians. 

Table 11-2 shows the uncertainty in lamp calibration as a function of 
wavelength in the range of interest, also supplied by Optronics Laboratories. 
This is the "reported uncertainty in the NBS 1973 irradiance scale along with the 
uncertainty associated with the transfer from the NBS stamdard to the issued 
200 Watt standard". There is another 1% uncertainty quoted by Optronics in the 
stability of the lamp. 

Because of the large (40 nm) half-height bandwidths of the filters used with 
this pyrometer, it is necessary to calculate the effective central wavelength of 
each filter, which is dependent on the spectral dependence of the other optical 
components, and on the temperature of the source. The other optical com- 
ponents have a weak spectral dependence relative tp the filters, and this does 
not change greatly from shot to shot as components are replaced. It is thus 
possible to use nominal values published in component specifications, and meas- 
ured values for the expendable components (the mirror and window) without 


Table II- 1 


Standard lamp uallbralion 


Wavelength (nra) 

Spectral Irradiamce^^^ 

400 

.438 

450 

.908 

500 

1.52 

550 

2.21 

600 

2.90 

650 

3.52 

700 

4.07 

750 

4.50 

800 

4.77 

900 

5.01 

1000 

4.92 


1. Spectral Irradianco in minro''ratts per (cm^ nanometer) ut a distance of 50 
cm. 






Table 11-2 

Spectral Irradianoe uiicertainties of standard lamp 


Wavelength 

(nm) 

NBS Uncertainty*^* 
(%) 

250 

2.5 

350 

1.4 

450 

1.2 

555 

1.0 

855 

0.75 

800 

0.9 

1300 

0.7 

1600 

0.7 


1. Reported uncertainty in the National Bureau of Standards 1973 irradiance 
scale on which the lamp calibration was based. 

2. Uncertainty associated with the transfer from the NBS standard, reported by 
Optronics Laboratories, Inc. 
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greatly aflecting the results. The wavelength- dependent parameters of the vari- 
ous optical components (Fig. 11-3) are: 

Rh(X) = reflectance of the aiuminized turning mirror; a typical mirror was meas- 
ured at the Lawrence Livermore standards laboratory. 

Tw(X) = transmittance of a typical glass port window; measured at Caltech with a 
Cary Model 17 spectrophotometer. 

Ti,(X) = transmittance of one Jeagers 35B1543 objective lens; measured at Cal- 
tech. 

Tbj(X) = transmittance of Ealing 22-8924 pellicle beamsplitter; from Ealing 
Optics catalog. 

Rgj^(X) = reflectance of Ealing 22-8924 pellicle beamsplitter; from Ealing Optics 
catalog 

Tb-(X) = transmittance of Melles-Griot 03BTF007 plate-type beamsplitter; from 

m 

Melles-Griot Optics catalog. 

Rb-(X) = reflectance of Melles-Griot 03BTF007 plate-type beamsplitter; from 
Melles-Griot Optics catalog. 

Tj(X) = transmittance of filter i; measured at Caltech. 

T](X) = transmittance of Melles-Griot 01LAG003 aspheric glass focusing lens; 
measured at Caltech 

Rj(X) = ^j(X)TJ(^) net spectral responsivity of channel i. 
r|(X) = net spectral reduction factor of channel i 
r,(A) = R«T»nTB,V,Rd. 
rzM = R«T»T=LRB,HBjT,Ri. 
ra(X) = R»T,I«i,Tb,Rb,T,Rb. 

rM = Hi|TwI«iJ!B,TB,TiRd- 


L 

Rd 

T| 

I 

R 

I 

Fig. II'S. Schematic of optical pyromc 
explained in text. 
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The mean wavelength of channel i Ls therefore 

m 

/r,(\)T,(\)XdX 

X,= -^ (U-1) 

/r,(X)T,(X)d\ 

0 

There vrill also be a shift in effective wavelength as a function of tempera- 
ture of the source, euid this should be included. To calculate the shift, we must 
assume a priori that the source radiates a Planck function. This assumption will 
also need to be made in order to reduce the data, so nothing is lest in making it 
now. The data will confirm the validity of the assumption, in the cases studied 
(Appendix 111). 

The mean effective wavelength as a function of temperature is 

•• 

/ r,(X)T,(X)r(T.X)XdX 

X, = -L , (II-2) 

/r,(X)T|(X)f(T,X)dX 

0 

where f(T,X) is the Planck function 

he 

f(T,\) a X-S(e - 1)“^ (II-3) 

Ihe mean effective wavelengths of the channel-filter combinations which 
have been used are plotted in Figs. 11-3 and II-4, The calibration lamp tempera- 
ture is about 2700'’K, so the effective wavelengths for calibration are those at 
2700°K. The appropriate spectral irradiance is found by interpolating between 
values given in Table 11-1. 
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Appendixlll 

Measured spectral radiances and calculated 
temperatures 




The measured spectral radiances used to calculate shock temperatures are 
reported here. In Tables III-l and III-2, spectral radiances obtained in seven 
shots conducted with the narrow-band six-channel optical pyrometer at 
Lawrence Livermore National Laboratory are tabulated. The design, calibration, 
and execution of experiments with this instrument are described by LyTzenga et, 
al. (1979) and Lyzenga (1980). Because narrow bandwidth (10 nm wide) Alters 
have been used with this instrument, the effective wavelengths are not a func- 
tion of temperature. Table III-3 gives the spectral radiances measured at Cal- 
tech using the wide-band four-channel optical pyrometer for Ave experiments. 
This table includes the effective wavelength of each channel, which is 
temperature-, and therefore shot-dependent. 

The temperature is calculated in all cases assuming a Planck distribution of 
radiated light 

fig 

where is the spectral radiance in WSr“^m~^ e is the emissivity, \ is 
wavelength in m. and T is absolute temperature in ®K. The constants are 
Cl = 1.191 X 10“*® and C 2 - 1.439 x 10“*. The measured spectral radiances and 
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the best fitting blackbody Planck curves (with eal.O) are plotted together in 
Figs. III-l through IIM2. 












Channel 


Table Iil-2 

Anorlhile glas» speulral radiance data 


AnST 


AnlT 


Nf 

450.2 

11.4 

±0.6 

507.9 

13.3 

±0.5 

545.1 

14.9 

±0.8 

598.0 

14.3 

±0.4 

850.0 

14.5 

±0.7 

792.0 

11.5 

±0.3 


450.2 


507.9 


545.1 


598.0 


650.0 


792.0 


1.92 

±0.15 

2.73 

±0.29 

3.24 

±0.13 

3.72 

±0.23 

4,09 

± 0.21 

4.34 

±0.32 


An3T 


1. Wavelength (X) in nm. 


1 . Spectral radiance (N^) in 10*® Wm“®Sr 













Table 111-3 


Anorlhile glass spectral radiance data 



An5T 

An6T 

AnTT 

An9T 

AnlOT 











v/iiGuiri"L 











x(») 

nP 

X 

Na 

X Na 

X 

Na 

X 

h\ 

1 

454.4 

0.0455 

454.8 

0.0137 

off scale 

453.5 

2.18 

453.0 

23,5 



±0,0078 


±0.0076 


±0.09 


±1.2 

2 

558.8 

0.165 

557.2 

0.0595 

of! scale 

603.2 

3.98 

802.8 

19,0 



±0.014 


±0.0105 


±0. 14 


±0.7 

3 

647.6 

0,314 

±0.028 

847.6 

0.0128 

±0.019 

off scale 

747.5 

4.25 

±0.15 

747.3 

19.0 

±1.6 

4 

748.3 

0.412 

748.4 

0.0226 

903.3 0.623 

903.3 

3.77 

903.3 

11.9 



±0.038 


±0.024 

' 

±0.026 


±0.17 


±0.5 


1. Wavelength (X) in nm. 

2, Spectral radiance (N;^) in 10*® Wm~®Sr“^ 


O) o 









































Spectral Radiance (lO'^W sr‘'iTf^) 
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Appendix IV 

Particle Velocity Experiment Design 


The purpose of the particle velocity experiment is to provide Hugoniot, 
sound speed, and release data in a single experiment. The experimental design 
is similar to those of Grady et xd. (1074), Larson and Anderson (1979), and Kondo 
et al. (1980), and makes use of electromagnetic particle velocity gauges (Dremin 
and Shvevdov, 1974). Gauges are oriented in a steady, uniform magnetic field 
such that the active element of the gauge, the magnetic field lines, and the 
direction of motion are all mutually perpendicular (Fig. IV- 1), ? 4 nd the gauge 
leads are parallel to the magnetic field. An electromotive force is induced along 
the length of gauge element, proportional to the velocity of the gauge. For a 
gauge with effective length L in a magnetic field B, the potential measured 
across the gauge leads is 

V(L)=DLUp(L) (lV-1) 

where Up is the particle velocity of the material in which the gauge is embedded. 

Gauges were photoetched from 10.0 jj,m thick copper foil with a 12.5 p,m 
thick polyamide (Kapton) film backing. The active elements of the gauges were 
0.9 to 1.0 cm long (Fig. IV-2). The precise effective gauge length (used in equa- 
tion IV-1) is dependent on the geometry of the gauge due to edge effects from 
the leads and the finite gauge width. It is about equal to the center-to-c enter 
lead distance of .95 cm. 





jOr 
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Rock samples were constructed by sawing rectangular slabs about 3.2 cm 
by 4.5 cm and grinding them to a uniform thickness of about 1.5 mm. 
Archimedeam densities were determined for each slab individually. Four slabs 
were bound together with epoxy, with three gauges at the interfaces and one at 
the free surface. The mean thickness of the glue at each interface was 9 /zm, as 
determined by measuring the total thickness of the completed target and sub- 
tracting from it the thicknesses of the individual component slabs and gauges. 
All four sets of gauge leads extended from the sample in the same direction, 
along the long axis of the slabs, such that the gauges were all parallel and even. 
Coaxial cables were soldered to the gauge leads, and were connected to oscillo- 
scopes for the shot. A typical rock target is pictured in Fig. IV-3. 

Water targets were constructed similarly, with the slabs replaced by 1.5 
mni thick fiberglass rings. The polyamide membranes were tightly stretched 
over the rings, and the rings were bonded together with epoxy. They were 
movmted on a 3,0 mm thick polycarbonate (Lexan) driver plate, and each cell- 
like space between the sets of gauges was injected with distilled water. The dis- 
tance between the gauges was determined by focusing a microscope on each 
gauge and recording the difference in focal positions, and correcting for the 
index of refraction the water between the gauges. A water target is pictured in 
Fig. IV-4. 

In Fig. 4-1, the entire particle velocity experimental assembly is shown 
schematically, with a rock target mounted on the center of the axis of a set of 
Helmholtz coils, which supply the magnetic field. To ensure mutual perpendicu- 
larity between the gauge, field, and particle velocity, the target is aligned by 
shining an alignment laser' beam down the gun barrel and adjusting the target 
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Fig. IV-3. Typical rock sample, with free-surface gauge visible. 
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Fig. lV-4. Typical water sample. 




until the beam is reflected back up the barrel by a mirror mounted on the tar- 
get front surface. 

The Helmholtz coils (Fig. IV-5) have a radius and separation of 14.1 cm, and 
are wound with four turns of gauge 10 copper wire. The magnetic field at the 
center point of the coil axis is constant to third order, and equal to 


B= 


8/i„lN 

( 5 ) 3 ^ 


(IV-2) 


where ( 1 q is the permeability constant, 1 is the current, N is the number of turns 
and r is the radius and distance between coils. The field is typically about 1.8 
kG. The current is supplied by a bank of 5 15 capacitors charged to 5.0 kV, 
and reaches its peak value of about 6.7 kA in about 60 /.ts (Fig. IV-6). Because 
the time scale of the experiment— which is determined by the shock and rarefac- 
tion transit times through the sample (w 2/i,s)— is short compared to the period 
of the capacitor-coil circuit (f« 250/zs), the timing can be controlled so that 
impact occurs at the peak current, when the field is eflectively constant. 

A simplified circuit diagram is shown in Fig lV-7. As the projectile 
approaches the end of the gun barrel, it interrupts a timing laser beam. This 
sends a signal to a trigger delay unit, which, after a pre-determined delay, sends 
a pulse to an EGScG TM-11 trigger module, closing a spark gap which in turn fires 
an ignitron, which dumps the charge in the large capacitor bank through the 
Helmholtz coils. The current is monitored by a Pearson Electronics Inc, Model 
301x "Viride Band current transformer, which is fed into a Tektronix Model 545A 
Oscilloscope, which is triggered by the same pulse that triggers the ignitron. 




Current (kA) 
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Fig. lV-6. Oscillogram showing time dependence of current through Helmholtz 
coils during experiment. Sharp pulse at current peak is fiducial mark, which 
indicates time at which impact occurs. 
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CURRENT CIRCUIT I SIGNAL CIRCUIT 


Fig. IV-7. Circuit diagram for particle velocity experiment. 
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In the signal circuit, self-shorting coaxial shorting pins are used to trigger 
all four oscilloscopes and a fiducial generator. The fiducial pulse arrives at each 
oscilloscope iOOns after the oscilloscope trigger. It is 5Uns wide, and is fed into 
the Z channel of each scope. It leaves a gap in the trace which allows synchroni- 
zation of the signals along the time axis. Each foil lead is connected to the 
center lead of a coaxial cable and the cable shields are grounded together. The 
two cables are connected to the two oscilloscope channels and terminated by 
son, and the scope is operated on differential mode in order to reduce line 
noise. Tekronix 475 and 485 models were used in conjunction with Polaroid cam- 
eras to record the signals. Time and voltage calibrations were conducted using 
the Tektronix 485 built-in calibrators. 

Not all the particle velocity experiments were conducted with the exact cir- 
cuit design and set-up outlined above. For example, in some experiments only 
one cable per gauge was used and the oscilloscopes were not operated in 
differential mode; and in the early (not always successful) experiments, smaller 
coils and different gauge fabrication techniques were used. The experiment 
described here was a product of evolution, and was the most successful. 
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